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Abstract 

Numerical simulation results (with emphasis on channel temperature) for a single gate AlGaN/GaN 
High Electron Mobility Transistor (HEMT) with either a sapphire or SiC substrate are presented. The 
static I — V characteristics, with concomitant channel temperatures ( T ch ) are calculated using the software 
package ATLAS, from Silvaco, Inc. An in-depth study of analytical (and previous numerical) methods for 
the determination of T ch in both single and multiple gate devices is also included. We develop a method 
for calculating T ch for the single gate device with the temperature dependence of the thermal conductivity 
of all material layers included. We also present a new method for determining the temperature on each 
gate in a multi-gate array. These models are compared with experimental results, and show good 
agreement. We demonstrate that one may obtain the channel temperature within an accuracy of ±10 °C in 
some cases. Comparisons between different approaches are given to show the limits, sensitivities, and 
needed approximations, for reasonable agreement with measurements. 

I. Introduction 

AlGaN/GaN HEMT amplifiers have achieved record power densities at microwave frequencies 
(refs. 1 and 2). Many of these “hero” devices have been small, and the aim of the experiments was to 
determine the power generation capability. The small units often consisted of one or two gates, so the 
problems of heat dissipation were not addressed. The large power density generated in power amplifiers 
causes considerable self-heating, and an accurate estimate for the channel temperature is needed for both 
design purposes and reliability estimates. For a practical device that will develop power at the multiple 
watt level, self-heating in the channel will be extreme. For power densities of 10 W/mm, the thermal 
management schemes will need to be more robust than those for other existing technologies, since the 
power density dealt with in AlGaN/GaN HEMTs is an order of magnitude larger. Knowledge of the 
strength of the self-heating is essential as increasing temperature reduces the mobility (it decreases as 
T 2 3 ), and other basic material parameters, with resulting decrease in DC and RF performance (ref. 3). 
Thus severe reduction occurs in the output power, transconductance, cutoff and maximum oscillation 
frequencies, and the reliability. Preliminary results (refs. 4 to 6) for device failure or permanent damage, 
have indicated failures occur when the channel temperature rises into the 200 to 300 °C range. However, 
reports of undamaged devices after operation at 600 °C have been reported (ref. 7). Both the long and 
short period output power changes measured in many devices under stress testing are related to the device 
temperature. The trade-off in the number of fingers in a MMIC cell is controlled to some extent by the 
thermal hot spots that develop near the edge of the gate in the gate-drain access regions in the center of a 
multi-finger cell. These factors show the need for reasonably accurate estimates of the channel 
temperature in multi-gate devices. 

The calculation of the channel temperature is, however, a difficult task for the following reasons. In 
principle, one should use a 2-D electro-thermal solver which performs a self-consistent calculation for 
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currents and temperatures for at least the basic cell (one source-gate-drain section). However, due to the 
large scale differences in the dimensions of the two-dimensional electron gas (2DEG) and the overall 
structure (substrate on the order of 350 pm), the simulations require complex meshing, and subsequently 
run for up to 15 hr. The 2-D self-consistent solver provides the initial estimate of the temperature field 
and the heat generation function g(x,y) at the chosen bias point. The 2-D slice is 1 pm thick and is 
assumed to be in the center portion (with respect to z) of the device. Next one should use a 3-D thermal 
simulator with g(x,y,z ) as the source function, to obtain a more refined estimate for the temperature field. 
However, the prescription of generating g(x,y,z ) from g(x,y) depends on the particular device geometry, 
etc. This approach has been emphasized in (ref. 8), as large inaccuracies from the 2-D runs may be 
present. Apparently the heat flow from the ends of the gate fingers can significantly affect the temperature 
field near the center (by up to 43 percent). One must also accurately define the boundary conditions; so 
the details of the top surface (passivation and metal layers) are needed. This information is usually not 
known precisely. This is a problem, as the heat generation region is sensitive to changes in physical 
details, bias, and the baseplate temperature. In any situation one must include the change in thermal 
conductivity with temperature of all material layers, or large errors will quickly accumulate. 

Several 2-D and 3-D simulators (refs. 9 to 12) are available to determine temperature profiles. A 2-D 
electro-thermal package, MINIMOS-NT (ref. 13) considers the coupling between the current and heat 
transport equations, and thereby performs a self-consistent calculation. The SILVACO (ref. 14) software 
package, which we use, is also a self-consistent 2-D electro-thermal simulator. 

In this memo we first present analytical models for channel temperature determination, and compare 
their predictions with measurements and simulations found in the literature. The accuracy is generally 
within 10 percent of measurements, and in some cases much better. These analytical methods are for both 
single and multiple gate devices, with any epitaxial layer sequence. We can handle effects of package 
dimensions, chip separations, and include the dependence of the thermal conductivity on temperature for 
all materials. Then we present results from the 2-D self-consistent electro-thermal solver ATLAS from 
SILVACO. We simulate I-V curves and measured temperatures from the only set of publications in the 
literature that have presented simultaneous I-V characteristics, and concomitant thermal measurements 
(refs. 15 to 18). We find good to excellent agreement. We have, for the first time, replicated measured 
I-V characteristics and channel temperatures for AlGaN/GaN HEMTs on both sapphire and SiC 
substrates. 

As stated earlier, the calculation of the channel temperature is not easy, and its measurement is also 
not easy. Achieving both accuracy and spatial resolution is difficult, so a reasonably simple analytical 
prediction for amplifier/device level design purposes can be very helpful. When using such 
approximations, it is helpful to know the sensitivity of the final temperature with the many parameters 
involved. This sensitivity is investigated here. A final comment concerns the uncertainty in many of the 
critical device parameters. One can question the use of sophisticated software packages, when some of the 
material parameters are unknown by as much as an order of magnitude. Unfortunately this is the case, as 
the critical parameter (thermal conductivity) depends on a given layer’s doping, growth procedure, 
process steps, and physical condition (cracks, etc.). This bolsters the need for a reasonably accurate 
analytical treatment. 


II. Previous Analytical Methods 

At the present time many investigators have studied the thermal properties of AlGaN/GaN HEMTs 
and GaN FETs under various conditions (refs. 19 to 31). They have presented combinations of measured 
data, 2-D and 3-D simulations, Monte Carlo results, and analytical schemes. The aim was to determine 
the channel temperature T ch as the dissipated power changed. The change in T ch may be due to changes in 
the bias point, baseplate temperature, or power added efficiency (PAE). In many cases constant thermal 
conductivity values were assumed, as this simplifies the problem (treats the system as linear). Figure 1 
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displays two possible plots of the channel temperature versus dissipated power. By definition, A T ch (the 
increase in channel temperature) and the dissipated power P diss are related by the “thermal resistance” R,/„ 

A T c h=Rth ' Pdiss (1) 


The final temperature T ch , is A T ch + reference temperature. The reference temperature is generally the 
baseplate’s value. From inspection of the figure one sees that R th is not constant, but varies with P diss . This 
variation is due to many factors. One is the change in thermal conductivity as the temperature increases (it 
decreases with increasing temperature), and another is the change in the heat generation function g(x,y,z ) 
as the dissipated power increases. Generally g(x,y,z ) increases in magnitude and the volume over which 
this region extends gets more sharply defined. Also, the location of the peaks of g(x,y,z ) change as the 
current path through the device changes with bias and RF drive. The nonlinearity may be concave up or 
down, depending on circumstances. These differences may be due to the geometry of the finger layout, 
the edge effects, or the rise in the baseplate temperature with increasing power dissipation. For example, 
for a sapphire substrate, the channel temperature may be in the 320 to 370 °C range for dissipation near 6 
W/ mm . This temperature may be about 1 00° higher than the extrapolation from the low temperature 
linear assumption. For SiC at this same dissipation, the temperature may be only 95 °C. But in both cases, 
one cannot extrapolate from the low temperature line. In general then, R th is a function that depends on 
many parameters, such as thermal conductivity and the detailed device geometry. This functional 
dependence will be explored in detail in subsequent sections. 

The determination of channel temperature from analytical expressions is very difficult at the present 
time for many reasons. One complication is the appropriate form for the thermal conductivities k(7) as a 
function of the temperature T, for the various material layers. While assuming constant k simplifies the 
analysis (treats the system as linear), we will show that the error is quite large for AlGaN/GaN HEMTs on 
either sapphire or SiC substrates. Even when this dependency is included, uncertainty exists for the proper 
values of k(7) in each layer. Since different growth methods and doping levels produce very different 
material properties, the reported measured values for k(7) have a large dispersion. For GaN (refs. 32 to 
37), for sapphire (refs. 38 to 40), and for SiC (refs. 41 to 44). A second problem is the choice of the 


Representative device temperatures versus dissipated power 



—♦—SiC substrate 

—■—sapphire substrate 

—♦—linear interpolation 
from 27 °C 


Dissipated power (W/mm) 
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Figure 1. — Representative channel temperatures versus dissipated power 
for devices on sapphire or SiC substrates. The secants to the curves 
(originating from T ch = 27 °C) are representative of the thermal 
resistances. The fact that the slopes are not constant reflects the 
changes in thermal resistance with dissipated power. 
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best approximate equations to use, as many have been published, and in many cases they don’t all 
converge to a common value. We have also found that the discrepancies are due in part to errors in 
published equations and mislabeling of plots. We list a few problem areas: 

1. In self-consistent 2-D simulators, does it matter that one uses the simpler drift-diffusion (D-D) 
approximation, or the more basic hydro-dynamic model? However, for the latter, what are the appropriate 
relaxation times? In the D-D formulation, one does not treat hot carrier effects (due to large electric 
fields) adequately. Also, what is the proper form for the mobility as a function of field, temperature, and 
doping? What is the effect of changing scatter-limited velocity with temperature? We assume some of 
these effects are minor, as will be shown later when we show excellent agreement with experiment. 

2. Shorter gates enhance self-heating, and this is a very strong function of gate length. 

3. Some experiments have shown variation of temperature along the gate fingers to be as much as 
100 °C, 

While many papers have appeared with methods to calculate the channel temperature, we will only 
concentrate on a few in this section. Recently Darwish, Bayba, and Hung (refs. 45 to 47) presented 
original closed-form expressions for calculating the thermal resistance of multi-finger AlGaN/GaN 
HEMTs and FETs. The two papers are complementary (one primarily for HEMTs, the other for FETs) in 
that somewhat different assumptions were used in each. The model consists of a vertical slice through a 
multi-fingered device, wherein the slice width is determined by inclusion of the source and drain pads 
with the gate in the center. The HEMT model specifies the GaN layer as region I, and the substrate is 
divided into two layers; regions II and III. In the first paper the isotherms were assumed to be circular 
cylinders, prolate spheroids, and elliptic cylinders, in regions I, II, and III, respectively. In the second 
paper, only a single layer is considered (appropriate for FETs), and the isotherms were assumed to be 
confocal ellipsoids and elliptical cylinders, as one moves from the gate to the heat sink. Attempts to verify 
the model included comparison with numerical simulation using ANSYS (for shapes of the isotherms and 
thermal resistance), and experimental data. The calculation assumes constant thermal conductivity and is 
restricted to multiple finger (greater than about 5) devices. They state the effects of non-constant thermal 
conductivity may be included by using the Kirchhoff transform; however they do not demonstrate the 
calculation. The wafer is assumed to be large in lateral extent, so end effects are neglected. The heat is 
generated by the gate electrode, and the thermal resistance is calculated assuming adiabatic boundary 
conditions on the vertical walls of the slice. This assumption implies that a few identical cells must exist 
on either side, so good results should be expected for devices with at least 5 or more gates. The 
demonstration of the close fit between the analytical equations and the shapes of the isotherms and 
thermal resistances as determined by the simulator ANSYS, makes the use of their results very 
compelling. There are, however, some drawbacks. The dependence on gate length is questionable; by this 
we mean the model for the gate shape is not completely convincing. The fact that the analytical equations 
results agree with the numerical simulations, means that the analytical model is correct. It does not mean 
that the model correctly reflects the heat source shape in actual HEMTs. The neglect of k(T) may have 
some effect on the shape of the isotherms as the temperature increases, and the resultant thermal 
resistance. In other words, the change in R th with dissipation level (due to changes in g(x,y,z) and k(T)) is 
omitted. Also the variation of R th with the number of fingers is eliminated by the nature of the calculation. 
Their equations are given below; see figure 2 for clarification of parameters. Note that Darwish et al. uses 
the symbol 0 for the thermal resistance R,i,. 

Rth (total) = R th (i) + R th (il) + R, h (ill) 
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Figure 2. — Geometry for calculating thermal resistance, 
as given by Darwish, Bayba, and Hung (refs. 45 to 47). 
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III. Thermal Resistance Formulas for Single and Multigate HEMTs 

In this memo we consider both single and multiple gate devices and include the dependence on 
temperature of the thermal conductivity as well as the actual size of the MMIC chip ( end effects). The 
general layer structure for a typical device is shown in figures 3 and 4. Figure 3 gives a top view of an 
eight gate device, while figure 4 shows the approximate layer thicknesses and thermal conductivity ranges 
for all materials. Our procedure to determine the temperature of the gate electrode (which we assume is 
nearly the same as the temperature in the GaN channel) is a series of about nine steps. In general we 
determine the thermal resistance of each layer of the device, and from that along with the I-V product of 
the bias point, the temperature is determined. The schematic in figure 5 gives the geometry for both single 
and multiple gate situations for the calculation of the thermal resistance, R,/,. The following list of steps 
gives a detailed explanation of our calculation procedure. 



Cu-tungsten carrier plate 


Figure 3. — Schematic of an eight gate AIGaN/GaN 
HEMT for a microwave power amplifier. 


0.1 -.175 pm 

f ' 3 N4 *■ i i k = 0.0096 

(passivation layer) 

0.02 pm 

GaN (cap layer) K = 1 .6 1 1 

0.022 pm 

AI0.2 Gao. 8 N " l 

(barrier layer) K | | 

0.5 !im 

GaN / n-, i6 1 I 

(buffer layer) / a i culated I | 

2 deg* — 7 1 j 

0.015 (im 

AIN (nucleation „ „ ■ ■ 

. \ K = 0.1 SI 

layer) s « 

100 pm 

Au source " " 

via k = 3.18 — 1 

SiC K fTo) = 3.4 1 1 

k(T) = calculated | _ _| 
35 pmi- h 

20 pm 

AuSn k = 2.4 

AuGe k = 0.9 

Epoxy k = 0.3-.6 

Silver epoxy k = 0.016-.075 

: 1000 pm 

Cu-tungsten 

(carrier) k = 4 


Note: The units for thermal conductivity k are W/(cm-K) 
* Two dimensional electron gas 


Figure 4. — Cross-section model of the FIEMT device, 
showing approximate layers of material with average, 
or ranges, of thermal conductivities. 
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(b) 


Figure 5. — Geometry for calculating thermal 
resistance for (a) a device with a single gate 
finger, (b) multi-gate structures. 


Thermal Analysis Procedure 

1. For a single gate device, start at the bottom surface and assume the baseplate temperature T h is 
known. Assume the temperature drop across the adhesive (solder) layer is small, so the bottom of the 
substrate is near T h . 

2. Calculate k (T h ) for either the sapphire or SiC substrate. 

3. Calculate R, h for the substrate. We recommend the formulas published by Masana (ref. 48), which 
are valid for a single heat source with the overall size of the chip considered. See figure 5(a) for the 
geometry. The equations are: 


Lg 21 x 

W„ = WtL x 
y<< 

Ys 

Acm / \ Av 

P, = K,/k,+i 
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vv n 1 


tan a = 11 - 


( " i+i 


P.9 I 


+ /„ 

1+p, n 


A 


tanP = 


V 


l-l ^ 
Y 


\ w„ +t^— ^y* 


i+p, n 




l , 

% + * x Ye 

1 + p, 11 


,/f ^ ^ 4 k(t)/ x y e tan a - tan P 


( 3 ) 


l x +w tan a 

, w tan P 
l x + 


4. In the previous step one needs the dimensions (2 l x , 2 l v ) of the effective flat strip heat source on 
the top surface of the substrate. The critical parameter is the effective gate length 2 l x . For SiC use 
4.5-(gate length), and for sapphire use 6.5-(gate length). There is some question as to the sensitivity of this 
assumption, as the actual size of the heat source is dependent on the bias point and the geometry of the 
structure. This is a critical step, as the final temperature is sensitive to the size of the heat source. All 
simulators that are not coupled electro-thermal calculators suffer with this problem, as they require the 
user to define the heat generation region g(x,y,z). The formula for R, h is for constant k, and the value is 
that for the temperature of the lower surface. With R th determined, calculate the linearized or “equivalent 
temperature” 0 at the top of the substrate (called the interface, int) 

Qint — Rth ' Rdiss 

where is the power dissipated for the single gate device, Pj iss = V d ■ I d , i.e., the product of the 
quiescent drain voltage and current. 

5. Find the actual temperature at the interface using the inverse Kirchhoff transform (discussed in 
more detail later) 


7int=^(eint) 

6. Calculate K(r int ) for the bottom of the GaN layer. 

7. Calculate R,/, for the GaN layer similar to the above, but now use 21 x = 1 . 125-(gate length) for 
sapphire or (gate length + gate-drain access) for SiC. 

8. Calculate the linearized temperature at the top of the GaN layer, which we assume (define) to be 
the channel value, 


°top = Rth (GaN ) • P ii ss 


9. Calculate the channel temperature 

r = ^ N (e, op ) 

where the inverse transform used is that which is appropriate for the form of the thermal conductivity 
variation in GaN. 

Observe p. s = k,(7)/ k,+i(7), for the i-th layer above the z-th +1 layer, thus the spreading angles a and P 
depend on the temperature T. We neglect the AIN layer which exists at the interface, as only a few 
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degrees exist across it in most cases. The effect of the low thermal conductivity of the AIN layer is 
reflected in the effective gate lengths used in the calculation. The low k of the AIN tends to funnel the 
heat entering the top surface from the higher conductivity GaN layer into a narrower region than a 
spherical spread would be expected to do, before it enters the substrate. 

For the multi-gate problem, we suggest the following procedure; first use the equations of either 
Darwish et al. (ref. 45), or Cooke (refs. 49 to 51, and 26). In these equation sets, use either the actual gate 
lengths for the heat source regions, or those indicated above. Next, use our more detailed calculation 
method; which will be described in detail in the following section. This procedure gives the temperature 
on each finger, unlike that of the above schemes, which can only give an “average” channel temperature. 
The formulas from (ref. 49), with typos corrected, are as follows (again refer to fig. 5): 


u = - 


M - 


cosh 

71 

(5 + 2O 


4 

l w 7 


cosh 

71 

f r\ 1 \ 

s-2l x 


4 

l w J 



2yju + 1 
yfu ~ 1 


P = 2 


cosh 

Tt 2^ 

4 w 

+ 1 

cosh 

k 2L 


X 

4 w 

-1 


Rth(T)= 


2/ v k(t)h: 

2(11 - 1) (77-2) 

In M InP 


( 4 ) 


We digress here to review the Kirchhoff transform. In principle, in any heat transport problem, one 
solves for the temperature field from the heat flow equation 

V-{k(t)VT}= -g(x,y,z), (5) 

where k( 7) is the temperature dependent thermal conductivity in the z-th layer, T is the temperature, and 
g(x,y,z) is the dissipated power density in W/cm 3 . This quantity is also called the heat generation density 
(source) region. The equation is altered by defining an equivalent, artificial, “linearized” temperature via 
the Kirchhoff transform (refs. 53 and 54), 


e = To+ W7)\ K ^ 


( 6 ) 


which relates the three temperatures T, T a , 0; the actual, reference, and linearized or transformed value, 
respectively. The transformed, linearized, or effective, temperature on the top surface of the substrate is 

02 =T b + R th (substrate) x P diss , (7) 
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where the known constant value for k (T 0 ) has been used in the determination of /G, (substrate). The real, or 
actual, or physical temperature at the top surface is determined from the inverse of the Kirchhoff 
transform 


7-2 =K S 7, 1(02 i (8) 

For the functional form 

K(r) = K(r o )(300/r)'\ (9) 

where T a = 300 K unless stated otherwise. The inverse Kirchhoff transform formula is 

^actual = {[6 - Ke - T 0 )\/T 0 r f (10) 

and the inverse is 

0 = l/(r-l)[rr o -r o 'Tl-'-]. (11) 

When r = 1 , we have 

^actual = T 0 exp([e -T 0 \/T 0 ), (12) 

where 0 is the effective temperature, and T a is the temperature at the bottom surface (the reference value). 
The temperature on the AIN top surface is 


T 3 =T 2 + Rfh (AIN) x (P diss /n), 


(13) 


where n is the number of gates. This assumes each gate dissipates the same power. Experimentally we 
know the gates near the middle are the hottest, and they probably dissipate the largest power. Finally, the 
junction, or channel temperature, at the AlGaN/GaN interface is found by 


0 4 = ^3 + Rfh (GaN) x {Pdiss /n) (14) 

^4 = ^GaN (®4 } • 05) 

Figure 6 summarizes the approximate heat flow pattern and the temperatures in the thermal resistance 
stack. Observe the AlGaN layer is neglected as it is extremely thin (up to maybe 500 A at most). 

At this stage of development of GaN technology, the actual values of the thermal conductivities of the 
various epitaxial layers, and their variations with temperature are not known to great accuracy. This 
necessitates that we must generate reasonable expressions to continue with the calculation. We choose the 
following expressions: 


GaN: 


Old 

k{t) 

New 

*{?) = 


300 

\ T ) 
300' 


(16) 

( 17 ) 
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Kuball’s (ref. 15) 


ic(r) = 1.6 


300 


(18) 




V J- J 


SiC. 

For doping levels near 10 17 cm 


For semi-insulating (S.I.) 



Sapphire : 


k (r) 


73.9 
T - 159 


0.49 


300 


(19) 

( 20 ) 


( 21 ) 


Figures 7 and 8 display the variations with temperature for GaN, and SiC, and sapphire, respectively. 
For GaN, the old expression was developed by assuming the coefficient was similar to that of GaP (which 
was measured), and averaging published room temperature values. The new equation utilizes some 
recently measured data. Observe at 300 Kthey differ by 18 percent. The expression by Kuball assumes a 
1 IT variation, and was used in his simulations. The value for SiC is an average of many experimental 
determinations. For sapphire, the first expression is a curve-fit of experimental data, while the second is 
an approximation for the basic 1 IT form. Figure 9 gives the variations of GaAs, Si, and InP; and the 
analytical expressions are: 


GaAs: 


k (T) = 0.47 


300 


xl-2 


( 22 ) 


Si: 


k(T ) = 1 .45 


( 300 xl-324 


: 1.45 


300 


A-2 


: 1.7 


f 30oV' 4 


(23) 


InP: 


k(t) = .658 


300 


\ 1.48 


(24) 


It is obvious that the value for all materials at room temperature, is uncertain by 8 to 18 percent, and 
the variation with temperature has some dispersion. The reasons are due to the quality of the tested films, 
and the range over which the fit was obtained. Thus we must be prepared to predict temperatures only to 
within 5 to 10 percent of measured values due to these uncertainties. 
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Figure 6. — Thermal resistance stack and schematic 
(not to scale) of heat flow through the layers. 


Thermal conductivity approximations for GaN 
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Figure 7. — Approximations for the temperature 

dependence of the thermal conductivity of GaN. 
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Thermal conductivity approximations for SiC and sapphire 



300 350 400 450 500 550 

Temperature (K) 

Figure 8. — Approximations for the thermal conductivities of SiC and sapphire. 
The upper-most curve is for semi-insulating SiC, while that below it is for 
doped SiC. The sapphire values are the lower two curves. 


Semiconductor thermal conductivities 



300 350 400 450 

Temperature (K) 

Figure 9. — Approximations for the thermal conductivities of Si, InP, and GaAs. 

Figure 10 gives an indication of the differences between 0 and T for our situation. The dotted curve is 
the actual temperature T versus 0 for the specific form for k(7) given on the plot. The form is that for SiC, 
the best heat conductor of all considered substrates. The reference temperature is T a = 300 K. A similar 
shaped (displaced) curve would occur for a different value for the reference. The lower line has unity 
slope that would correspond to the actual temperature equaling the transformed value (i.e., constant k). It 
is there to aid the eye to show how rapidly the actual temperature climbs as the equivalent one increases. 
The increase above the T = 0 line shows the effect of the decrease of k(T) from its value at T„. For 
example, the conductivity at 300 K is 3.3 W/(cm-K). If the conductivity were constant at this value, the 
thermal resistance would be calculated with this value for k. For a given heat flux, suppose one calculates 
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Actual versus equivalent temperature 
600, 1 



to 

€ 10 °- 
< 


0 - - - I 

300 320 340 360 380 400 420 440 460 

Equivalent temperature (K) 

Figure 10. — Actual ( 7) and equivalent (9) temperatures versus 0. The equivalent 
temperature graph is a straight line. The actual temperature curve (upper-most) 
is curved. The reference temperature is 300 K. 

9 to be 450 K, then the actual value would be 534 K (see the plot). Observe the increase is 84 K, which is 
extremely large! This shows the need to include the variation of k (T) in the calculation. One difficulty 
with the transform method is that jumps in the effective temperature occur at material interfaces, unless 
the forms for k( T) satisfy certain constraints (ref. 54). The constraint is that the ratio of conductivities of 
adjacent layers must be independent of temperature. In other words, the nonlinearity in the basic heat 
equation is removed with the transform, but makes itself felt, and must be dealt with, in the boundary 
conditions. For our case of multiple layers and sequential calculation steps, we calculate the actual 
temperature in each layer, and thus avoid any complications. 

The following table provides a compilation of reported values for thermal conductivities for 
components that are of interest to one studying the field of high temperature amplifiers. 

TABLE 1.— THERMAL CONDUCTIVITIES FOR 
MATERIALS USED IN HIGH POWER AMPLIFIERS 


Material 

K(W/(cmK)) 

GaAs 

k (7) = k (300)(300/7)‘- 2 k (300) = 0.44 
= 745 r 1 - 3 

= 56 r°' 87 3 50 K< J<500K 

SiC 

3.3(TJT ) l 5 (doped) 
~ 5 for S.I. 

Si 3 N 4 

0.016 

Sapphire (A1 2 0 3 ) 

k (300) = 0.25 - 0.5 
k (7) = k (300)(777),) 0 ' 8 

Au 

2.95 to 3.18 

AuSn solder 

0.573 

Silver epoxy 

0.0714 to 0.075 

Epoxy cements 

0.0188 to 0.0376 

Solders 

1.88 to 2.4 

Copper 

4 

Si 

K(r) = K(300)(77r o ) L324 
= 1.54(300/7) 4/3 
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Figure 1 1 . — Geometry for multiple gate (heat sources) 
on a substrate. 


IV. Multiple Gate Thermal Resistance Via Multiconductor Coupling Model 

The previous sections presented the accurate model from Masana for single stripes, and the two 
methods for multiple gates (Cooke and Darwish et al.). The equations of Cooke yield a single 
temperature, which is apparently an average value of all fingers involved. The formula is dependent on 
the number of gate fingers. Those of Darwish also give an average temperature, which is apparently that 
of the central region. The equation is assumed valid for a device with at least five fingers, as the number 
of fingers is not considered. 

Since the layout of gates is similar to the arrangement of multiple conductors above a ground plane, 
concepts from multiple conductor theory should be applicable for our problem. The problem to be solved 
is the determination of R t i, and temperature for each gate electrode as shown in figure 1 1 . In general, the 
Source-Drain (S-D) spacing is much smaller than the separation between each S-D region, so the heat 
generated in each GaN section can be assumed to not couple to heat generated by adjacent sections in the 
uppermost part of the device. Therefore, we are considering rectangular sources (the appropriate 
footprints) on the substrate top surface. Once the temperature on this surface is determined, then similar 
calculations, as given previously, will determine the temperature at the channel. Recall the equation for 
steady state heat conduction (assume constant k) is 

?2r= _ g(w) , (25) 

K 

Q J y Z \ 

Where g(x,y,z ) is the heat generation function. Poisson’s equation is V 2 V = - ’ " ’ , and the 

G 

similarities are obvious. 

We will use developments in electrical theory to ascertain corresponding results for temperature. This 
is the same motivation utilized in Cooke’s formula; his model uses the capacity for strips above a ground 
plane. For a finite number of conductors above a ground plane we know the charges and potentials are 
related by the Maxwell capacitance matrix [C]. 
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where the matrix [P] is the inverse of [C], and is called the potential coefficient matrix. If the conductors 
are assumed to propagate a pure TEM mode, then an inductance matrix is defined where 

[L][C] = » 0 e„[l] (26) 


or 


Wb]=Eo 




or 


[P] = — [L] (27) 

Po e 0 

Now we make the analogy V<r^T, where Pdi SS is the dissipated power (no relation to the 

potential coefficient matrix). Now assume py+^Ry, where Ry are the self and mutual coupling thermal 
resistances of the gates; then 


" T\ " 


1 

>3 

K> 

■ ■ Rm 

PrfLSSI 

t 2 

— 

Rn 


PclISS2 

/n_ 


_Rm 

■ ■ 7?nn. 

PdlSSN _ 


We assume all gates dissipate the same power, so the temperature on any gate is 

Tk = (Rn + Rk2 + Rm +• ■ ■+ Rktd Pdtssk (29) 

or just the sum of all thermal resistances in row k. Here Pdis.sk is the dissipated power in gate k, which is 
the same for all gates. We then assume 


_ Pu_ 
R\2 P\2 


(30) 
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that is, the ratios of the self and mutual terms are the same. This should be exact due to the analogy. We 
will determine Ry for any single conductor; then scale all other R,, accordingly. For our typical cases the 
parameter ranges are: 

W~ 24 pm /; ~ 350 |um Wlh~l- 10 -2 S/W-2 

S ~ 40 pm / ~ 200 (am l/h ~ 0.6 S/h ~ 0.1, 

which are not typical for most microstrip coupled line structures. 

Notice the thermal conductivity k and s correspond in the analogy. The analogy holds for 
homogenous materials, but for the inhomogeneous case, problems arise. Our case is inhomogeneous as 
we are assuming strips on a material with air above in the electrical case. The thermal case is that for 
finite thermal conductivity below the strips, but zero above. The electrical case cannot have zero for the 
permittivity of the air. Recall that microstrip is quasi-TEM and an effective permittivity s re is used. There 
are flux lines and potential contours in the air above the substrate, and & re is a measure of this effect. For 
the electrical case we have 


p 0 s 0 s re [p] = [L], (31) 

and thus we might replace k of the substrate with some effective value. We don’t have to deal with this 
however, as we are using ratios, and all terms carry the same k. Actually the gate footprints on the 
substrate are about 1 pm below the actual air (i.e., beneath the GaN and AIN layers), so we can assume 
we are dealing with an almost completely homogeneous region. Then the analogy will hold, i.e., we are 
solving Laplace’s equation in a source-free homogenous region. Recall that we use k( T„) to treat the 
problem as linear, then obtain the true temperature via the Kirchhoff transform. 

Our task is to determine the self and mutual inductance entries in the [L] matrix 



Ai 

L\ 2 ■ 


An 

M- 

L 2l 






_AsT1 




Amn _ 

and for say a 5><5 

array we 

have 




Ai 

A 2 

A3 

Lu 

As 


A 2 

A, 

A 2 

A 3 

Lu 

Fl= 

L u 

A 2 

Ai 

A 2 

A3 


L\a 

A 3 

A 2 

Ai 

A 2 


As 

Lu 

A3 

A 2 

Ai 


(32) 


(33) 


Observe all entries on any diagonal (top left to bottom right) are the same. The actual determination of the 
Lij is a complex problem. For one thing, inductance requires closed loops, whereas the [L] entries which 
are in units of Him are distributed values for assumed infinitely long lines. We choose two alternate 
methods to calculate \L\, the distributed, and the partial inductance approach. The partial inductances 
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depend on the length of the conductors as opposed to the distributed idea for infinitely long conductors. 
Therefore we are including the 3-D effect for heat emerging from the ends of the fingers. 

For N circular wires spaced S meters apart, the self inductance La and mutual inductances Ly 
(distributed) are (refs. 55 and 56) 


Lii=^-£n\ 
2n 


r 2h\ 


h))a , 


y u j 


where a = radius, h = height above ground plane, and 


m | Sfj + (lh ) 2 

1 4n I 5? 


For a 5 conductor array we have, S\ 2 = S, Sn = 25, 5| 4 = 35, Sis = 45; thus 

4n S 2 


(34) 


(35) 


(36) 


h3=—tn- 


Fo J( 2S f + ( 2h f 


[ (2 5) 2 J 

Due to symmetry of \L\, we need only generate the first row terms. Observe 

4n \ [(*-i>F ) 

For rectangular conductors: 


(37) 


(38) 


(39) 


(40) 


Lii ~lV n X + ln {w + t)\ 


where w = width of stripe, t = thickness, and 


L;: = -in- 
3 4rc 


5? + 




S ?J 


(41) 


( 42 ) 


We can also use the self inductance of microstrip to express La, 
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All of the above formulas are for external inductances, which should be appropriate for our case of 
widely separated gates. The previous formulas may be found scattered in the literature with slight 
modifications. One basic assumption is that Ly, the mutual inductance between conductors i and j, 
depends only on the distance between them. The presence of the other conductors is ignored. 
Experimentally this has been verified, and for our case of not strong coupling, it should hold well. 
Numerical solutions in the literature shows the errors incurred with this assumption (ref. 56); and for our 
parameter ranges the errors are a few percent at most. 

As stated earlier, our lines are not extremely long with respect to spacing and ground plane distance, 
so perhaps partial inductances could make for a better model. We choose only two forms; those for self 
and mutual inductances of rectangular conductors and filaments, respectively. 

L p u = 0.2 £ £n ( + 0.5 + .2235^-^ pH (44) 

l w + 1 ) i 



Then the entries in the \L\ distributed matrix are (ref. 57) 

i, j not equal to k. 


where the superscript p means partial inductance value. For example see figure 12. 

'{ L p n -L p n -L p l+ L p ,) (Zf 2 -Lf 3 -Tf 1+ Zf 3 )' 

M= 

_(^21 “ ^23 - ^31 + ^33 ) (^22 “ ^23 - ^32 + ^33 )_ 


(45) 


(46) 


(47) 


For the ground plane partial self inductance T 33 , we will use the image of a single conductor, i.e., 1 
or 2. 



Figure 12. — Notation for calculation of self and 
mutual partial inductances for two conductors 
over a ground plane. 
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All of the equations presented in this section have been checked with reports in the literature (refs. 58 
to 65), and we found good agreement in most cases. The worst cases were numerical simulations with 
tight coupling (not our case). For geometries close to our range (ref. 56), the error was ~3 percent. 


V. Comparison of Analytical Models — Single Finger 

This section compares the theoretical predictions for thermal resistance and channel temperature for 
single gate (finger) devices with the measurements of Kuball (ref. 1 5). 

One Stripe on Sapphire 

Kuball’s group simultaneously measured the channel temperature using Raman spectroscopy and the 
1 — V characteristics of AlGaN/GaN HEMTs. Figure 13 gives the schematic of a single 200 pm wide gate 
on a GaN layer grown on sapphire. The gate length is 4 pm. The source and drain pads are omitted for 
clarity, and their shape and presence did not have any affect on the channel temperature, as stated by the 
authors. The measured I-V characteristics are also shown. 

Table 2 summarizes the measured channel temperature, the drain current, and the dissipated power, 
for several drain voltages V d , all for the gate at 0 V. Notice that the dispersion in the data is from 10° to 
20°. At the bottom of figure 13 is a schematic of the gate region with the temperatures at several points 
determined by simulation (which will be discussed in detail in section VIII) at V d = 20 V. Observe the 
simulation values fall into the measured range of 179 to 190 °C. 


TABLE 2.— DATA FOR A SINGLE STRIPE HEMT ON A SAPPHIRE SUBSTRATE 


Er,V 

4 

8 

12 

20 

r, °c 

40 to 60 

80 to 95 

120 to 140 

179 to 190 

< T > avg. temperature 

45 

90 

129 

180 

I d , mA, drain current 

39 

42.5 

39.5 

32.5 

Pdiss, W, dissipated power 

0.156 

0.34 

0.474 

0.65 



Figure 13. — Schematic of single gate HEMT and the measured l-V curves. Notice the 
gate is not centered on the chip. 
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The procedure to calculate the channel temperature is simplified by noting that over 95 percent of the 
thermal drop occurs in the sapphire substrate due to its thickness. The sapphire constitutes 350/351.2 = 
99.7 percent of the thickness of the device. The key step is the determination of the effective gate length 
footprint on the top surface of the sapphire. For this case the actual gate length is 4 pm. If we assume the 
heat is generated uniformly under the gate, and it diffuses downward at 45° angles from the gate edges, 
the footprint at the bottom of the GaN layer is 6.4 pm. The assumption of 45° implies no edge effects, so 
the heat moves away into an unbounded semi-infinite half-space. Now passing through the AIN region 
causes some change in the heat flow, and we obtain best results by assuming the footprint on the sapphire 
to be 6 pm. This focusing effect in the AIN is demonstrated in section VIII. The parameters in section III, 
the Thermal Analysis Procedure, section II, for this case are: 

2 4 = 6 pm s T = 0.943 

21 y = 200 pm w = 350 pm 


2L X =180 pm 


2 L y = 300 pm. 


Then 


W„ = 3.89 


Ye = 33.3 
A =1.67 


l xn = 0.033 p s = 0.49/4 = 0.123, 


where we assume the solder layer on the copper heatsink is so thin as to be negligible. The thermal 
conductivity of copper is taken to be 4 W/(cm K). We find a = 43.8°, P = 15.7°, and R,h = 174.3 K/W. 
Table 3 provides the sensitivity of R t h for several effective footprints. 


TABLE 3.— -THE SENSITIVITY OF R th AND THE SPREADING ANGLES 
a AND |I WITH THE FOOTPRINT ON THE SAPPHIRE LAYER 


2/r (pm) 

0.5 

1.25 

2 

4 

6 

8 

R,h (K/W) 

294 

248.6 

226 

193 

174.3 

161.7 

an 

44.9 

44.8 

44.6 

44.24 

43.8 

43.5 

P(°) 

15 

15.1 

15.7 

15.7 

15.7 

15.7 


It will be shown later that only a few degrees are dropped across the GaN and AIN layers, so we can 
estimate R,u (sapphire) from the measured data. The procedure is the following. For the 4 V case the 
measured average channel temperature was 45 °C (318 K). Assume 2 °C dropped across the GaN and 
AIN layers, then the top of the substrate is at 3 1 6 K. Then the corresponding transformed temperature is 
therefore 315.6 K (from eq. (12)). The drop in the transformed temperature across the sapphire is thus 
315.6 K- 300 K = 15.6 K. Then R th (4 V) = AQ/P diss = 1 5.6 K/0. 1 56 W = 1 00 K/W. Table 4 gives the 
results for all drain values. 


TABLE 4.— THERMAL RESISTANCE AT SELECTED DRAIN VOLTAGES 


EXTRACTED FROM MEASURED RESULTS OF KUBALL 


Trf (V) 

4 

8 

12 

20 

R th (K/W) 

100 

163 

182 

188 
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Observe the average for the last three entries is 111 .1 which is slightly less than 2 percent different from 
our calculated value. Using R th = 174.3 K/W we find the channel temperature at V d = 20V as follows. The 
linear drop across the sapphire is A9 = (174.3) (0.65) = 1 13.3 K, so the top surface is at 300 + 1 13.3 K = 
413.3 K. The actual temperature from equation (12) is 437.7 K (164.7 °C). For the 1.2 pm GaN channel 
with 2 l x = 4 pm and using K(GaN) at 438 K, we obtain k = 1.97(300/438)° 48 = 1.64. Then the thermal 
resistance of the GaN layer is R,h = 6.1 K/W, and the change across the GaN is A0 = (6.1) (0.65) = 3.79 K. 
Using equation (10) with T a = 437.7 K, 0 = 441.67 K, and r = 0.48 yields 441.6 K, or 168.6 °C. So the 
drop across the GaN layer is just 3.9° as anticipated. This final temperature compares favorably with the 
average of 180 °C of the measurements. 

Table 5 summarizes the estimated temperatures for drain voltages of 4, 8, 12, and 20 V. Figure 14 is 
reproduced from Kuball’s paper and shows the temperature measurements in the access regions about the 
gate. The circles to the right of the data in the gate - drain access region are the values given in table 5. 
The asterisks are the values in the channel region obtained by simulation. The uncertainty in Kuball’s data 
was ±5 °C. We conclude that the outlined procedure works very well. 

TABLE 5.— COMPARISON BETWEEN AVERAGE MEASURED CHANNEL 
TEMPERATURES <T> AND CALCULATED ONES T calc FOR 


SINGLE STRIPE HEMT ON SAPPHIRE SUBSTRATE 



ov 


2q i i i i i i i i 

-2 0 2 4 6 8 10 12 

Position, pirn 

Figure 14. — The measured and calculated temperatures in 
the S-G and G-D access regions. The simulated values 
at the AIGaN/GaN junction are shown as asterisks at 4, 
8, 12, and 20 V. The calculated values are the circles. 
The uncertainty of any measured value is 10 °C. 
(sapphire substrate). 
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At this point we discuss some of the problems associated with this series of publications. For the 
results in these papers, the heat sink was not held at room temperature, but allowed to rise as the 
dissipated power increased. Thus the measured channel temperature value includes the thermal rises 
across the layers along with the u nk nown reference temperature of the sink. For the sapphire case the 
edge went up about 20 °C. Kuball states we can assume the sink rose by that amount. Our calculated 
temperature is about 12° low, which we can surmise is due to the fact that we assumed the sink to remain 
at room temperature (27 °C). The data in (ref. 18) suggest that the sink rose to 40 °C, which would boost 
our values approximately 13° ! Then we would agree to within 1°, which of course is beyond our 
anticipated accuracy. Nevertheless, we can apparently get rather good results with judicious 
approximations. Kuball states his system accuracy to be within 10 °C. 

Now k a1n is temperature independent and we assume its value to be 0.1. This is about 1/1 0 th of its 
bulk value, but typical for thin films grown on sapphire or SiC. As a matter of fact, the authors in (ref. 66) 
have reported measurements that show the thermal conductivity may be two orders of magnitude below 
the bulk value, when grown on Silicon. Kuzmik et al. (ref. 67) have used different values for k A in in 
simulations, depending on the layer upon which it was grown. Since we do not know the specifics of their 
AIN layer, some small error exists. However, the thermal resistance formula is insensitive to p.„ so the 
error should be very small. Assuming a parallel isotherm flow pattern in the AIN, (and the thickness to be 
0.015 pm) we obtain R,/, (AIN) = 0.75 K/W, with a corresponding thermal drop of 0.49°; which we have 
neglected. 

While we have been able to get results that agree with Kuball’s data, we still have some reservations. 
The first issue is the reliability of the measurements. In a recent paper (ref. 68) concerning the 
measurement of temperature in AlGaN/GaN devices using Raman scattering, it was shown that very 
different results are obtained when different phonon lines are observed. Specifically, the UV line recorded 
temperatures which were consistently higher by as much as 70 °C as that found using a line in the visible. 
The reason is the depth of penetration of the light into the sample. For a device similar to the one 
investigated here, the difference in measured temperature was about 70 °C for the dissipated power of 
0.65 W. Kuball’s group used the line in the visible, which averages the temperature in the GaN. 

Table 6 compares the thermal resistances calculated by the equations of Masana along with the 
corresponding channel temperatures. In these calculations we used the following approach. First start with 
the actual gate length and calculate the thermal resistance using the room temperature value for thermal 
conductivity. Use the calculated spreading angles to determine the footprint. Calculate the substrate 
thermal resistance. Then calculate the temperature rise across both layers using the dissipated power. Next 
recalculate R th using the values of k corresponding to the average of its values at the new temperature and 
room temperature. Then calculate the channel temperature. 


TABLE 6.— THERMAL RESISTANCES AND RESULTING CHANNEL 
TEMPERATURES USING THE EQUATIONS OF MASANA 


v* 

V 

P dissr* 

w 

Rth (GaN), 
K/W 

Rth (sapp), 
K/W 

T c h > 
°C 

4 

0.156 

5.37 

165.4 

53.6 

8 

0.34 

5.62 

172.3 

87.5 

12 

0.474 

5.81 

176.8 

113.6 

20 

0.65 

6.06 

182.6 

149.6 


Table 7 shows the results of a similar calculation, this time using the equations of Cooke. In this case 
we assume the spreading angles are 45°, since the edges are sufficiently far away. The final footprints in 
the two calculations are 4.064x203.15 versus 6.4x202.4 pm 2 . The final channel temperatures are 149.6 
and 93.4 °C, a difference of 56.2°, which we consider significant. Observe they differ from that in table 5 
by unacceptable amounts. 
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TABLE 7.— THERMAL RESISTANCES AND RESULTING CHANNEL 
TEMPERATURES USING THE EQUATIONS OF COOKE 


V d , 

V 

P dissi 

w 

Rth (GaN), 
K/W 

Rth (sapp), 
K/W 

T 

1 chi 

°c 

4 

0.156 

3.47 

92.5 

42 

8 

0.34 

3.56 

94.6 

60.4 

12 

0.474 

3.63 

96.5 

74.8 

20 

0.65 

3.72 

98.4 

93.4 


VI. Comparison Between Darwish et al., and the GRC-1 Model 

The previous section compared the predictions of Masana and Cooke for the thermal resistance of a 
single gate. Masana’s equations are only valid for a single gate, whereas those for Cooke are valid for 
either single or multiple gates. The formulas by Darwish, et al. (ref. 45) are only valid for the multiple 
gate case, and should be accurate for about five or more gates. The multi-conductor method of section IV 
should be applicable to either single or multiple gate cases. For brevity we will call the analysis method 
and final equation by Darwish, the one-shot result. The iterative approach, which uses Masana’s 
equations for the single gate where appropriate (in the GaN layer), and Cooke’s equation for the multiple 
gate case (in the substrate) will be labeled the GRC-1 method. The results from the multi-conductor 
equations will just be called the multi-conductor results. 

As stated previously the models are complementary, but they can be compared in their regions of 
overlap. We will determine a value for R th for the following data set: 

L g = 2 l x = 0.25 pm 
Wg = 21 y = 250 pm 

?GaN = 1 pm 

t S ic = 100 pm 
S = 40 pm 
K S ic = 3.3 W/(cm-K) 

KGaN = 1.5 W/(cm-K). 

which is that given in (ref. 45) as an example, and plotted in their figure 4. The one-shot method finds the 
thermal resistance of the GaN layer (region 1) to be 07=13.82 K/W, and that for the two portions of region 2 
to be 19.37 and 19.95 K/W. The total resistance is 53.1 K/W, which agrees exactly with the plotted curve in 
their figure 4, which is reproduced here as figure 15(a). The figure shows nearly perfect correlation between 
the analytical calculation and the finite element results. The GRC- 1 method obtains R,/, in the GaN to be 
30 K/W. This is assumed to be valid as the gates are separated by 40 pm, and will not interact in the 1 pm 
GaN layer. We have some leeway on the calculation of the thermal resistance of the SiC layer. Since the 
thickness is 100 pm and the spacing between fingers is 40 pm, and if we assume spreading is at 45°, then 
each finger interacts with four of its neighbors (two on each side). Using Cooke’s formula, the resistance of 
the SiC layer changes from R lh = 21 to 28.7 K/W, as the number of fingers varies from 1 to 20. Then the total 
resistance is 30 +(21 thru 28.7) = 51 thru 58.7 K/W. If we use Masana’s equation for a single finger, we 
obtain R th =24 K/W. The total resistance is now 30 + 24 = 54 K/W, see figure 15(b), which agrees closely 
with the previously found value from Darwish (53.1 K/W). While we hope for agreement, this result is 
surprising, and we are not sure if this is just a coincidence for this particular case. The fact that they give 
different values for the GaN and SiC layers, yet the total is the same, is unclear. The determination of the 
channel temperature is different for the one-shot versus the other two methods. For example; the 
temperature change for the one-shot case is (R t h)(P<im)(n) where P diss is the power dissipated per gate, and n 
is the number of gates. For n varying from 5 to 15 with P diss assumed to be 1 .0 W, the channel temperature is 
shown in figure 15(c). The temperature varies from 286 to 817 °C. 
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Reproduced figure from reference 45 




///////// 


Total = 53.1 K/W 
(b) (one shot) 


Tint 


■30 


■21-28.7 


///////// 

Total = 51-58.7 
(GRC-1) 



Total = 54 
(Masana) 


Substrate = SiC 



Figure 15. — (a) Reproduced figure from reference 45; showing the excellent agreement 
between the finite element code and their model’s results, (b) Observe that all three 
models give the same net thermal resistance; however, the distribution over the 
GaN and SiC layers in the one-shot case is very different, (c) The channel 
temperature from the one-shot calculation. 
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The GRC-1 and Mansana (only) method calculates the channel temperatures as follows. Assume the 
dissipated power is 4 W/mm, and each finger is 0.25 mm wide; thus the power dissipated per finger is 
1.0 W. The temperature change across the channel A T ch = /t,/,(GaN)( 1 .0 W), and that across the substrate 
is R th (S\C)(P diss )(n). The final channel temperature T ch is the sum of the changes across the substrate and 
channel +20 °C. Figure 16 gives the interface temperature as the number of fingers varies from 1 to 16, 
while figure 1 7 gives the corresponding channel temperature. The upper curve, (one-shot), is derived 
from the one-shot calculation, while the lower two curves are for Masana’s equations in the GaN and 
Cooke’s in the SiC, or Masana’s in both layers (neglecting any crosstalk in the SiC layer) respectively. 
Observe that this channel temperature calculation procedure gives the range of 270 to 663 °C as the 
number of fingers varies from 6 to 16 for the one-shot case. The upper temperature is around 150 °C 
lower for this case than the previous method. Table 8 summarizes the temperatures at the GaN/SiC 
interface and the channel, which are shown in figures 16 and 17. 



- Darwish 

- Masana 
GRC-1 


Figure 16. — Interface temperature rise as the number of fingers increases. The base 
is held at 20 °C. 



—♦-One-shot 
— ■ — GRC-1 
— Masana 


Figure 17. — Corresponding increase in channel temperature with change in the 
number of fingers. 
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TABLE 8.— INTERFACE AND CHANNEL TEMPERATURES (°C) FOR A GIVEN 
DATA SET AS COMPUTED BY THREE DIFFERENT METHODS. THEY ARE: 
MAS ANA ONLY (M), GRC-1, AND ONE-SHOT RESPECTIVELY 


n 

T'mi 

Tnt 

T\nt 

Tch 

Tch 

T c h 


M 

GRC-1 

One-shot 

M 

GRC-1 

One-shot 

i 

44 

41 

59 

74 

71 

73 

2 

68 

69 

99 

98 

99 

112 

4 

116 

127 

111 

146 

157 

191 

5 

140 

156 

216 

170 

186 

230 

6 

164 

185 

256 

194 

215 

270 

8 

212 

243 

334 

242 

273 

348 

10 

260 

302 

413 

290 

332 

427 

12 

308 

361 

492 

338 

391 

505 

16 

404 

478 

649 

434 

508 

663 

20 

500 

594 

806 

530 

624 

820 


Before considering more comparisons with published results, we discuss the various functional 
dependencies upon which the thermal resistance depends. One may write the thermal resistance as 
follows. 


R,k =Ato aN, 4„b, K GaN (70, k su \T), w g , l„ n, S) 

The eight parameters are the GaN and substrate thicknesses, their thermal conductivities, the gate width 
(per finger), the gate length, the number of fingers, and their separation. In general, R th decreases with 
increasing W g , increasing S, and increasing l x . It increases with increasing substrate thickness, and with 
the number of fingers, n. The strength of the variations depends on the geometry and the heat generation 
function g(x,y,z). 

The paper by Bertilsson et al. (ref. 8) concerning lattice heating in SiC MESFETs is illuminating for 
our purposes. Their figure 4 gives normalized thermal resistance versus finger width W g , n, and S. The 
reference values are 200, 10, and 100 pm respectively. Figure 18 gives R, h (normalized to the gate width) 
as a function of the number of fingers, with the gate length (2 l x ) as a parameter. The physical gate length 
was 0.5 pm, and the active layers constituted 1.38 pm. Using the footprint method, we determine the 
effective gate length on the SiC substrate to be 3.26 pm. The curve labeled reference 8 is extracted 
from figure 4 of the paper. The curve below it is that derived using the GRC- 1 method, with gate length 
= 3.26 pm. We consider the agreement to be excellent. The two upper curves show the result if the actual 
gate length, or the gate-drain spacing is used. This reveals the sensitivity to gate length as well as the 
agreement obtained using the actual length and determining its footprint. For the given gate length of 
0.5 pm, the one-shot methods gives the value of 9.6 K-mm/W. The agreement of both methods with the 
reference 8 results is aided in part by the thermal conductivity being the same in all regions; the 1.38 pm 
epitaxial layer and the 330 pm substrate. Figure 19 gives the variation of the thermal resistance with gate 
pitch 5; the diamonds are from reference 8, and the squares are the GRC-1 results. Some values obtained 
with the one-shot method are stated in the caption. Figure 20 depicts the dependence on gate width; the 
triangles show the GRC-1 result, and one-shot values are stated in the caption. The denormalized thermal 
resistance curve shows the true decrease with gate width. The one-shot results match the normalized 
curve closely, but the GRC- 1 curve is too flat. 
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Thermal resistance versus number of fingers 



— ♦ — 2/* = .5 
— ■ — 2/ x = 1 .2 
—a— Reference 8 
— x — 2/ x = 3.26 


Figure 1 8. — R* versus the number of gate fingers for different gate lengths (2 l x ), 
based on simulations by Bertilsson (ref. 8). The value predicted by Darwish is 
9.6 K-mm/W (independent of the number of fingers). 


Thermal resistance (normalized) n = 10, 2I X - 3.26 pm 



Gate pitch (S) (pm) 

Figure 19. — Variation of R th (normalized) with gate pitch. The values by Darwish are 
28.1, 9.6, and 7.7 at S = 20, 100, and 200 pm, respectively. 
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Normalized and denormalized thermal resistance versus gate width 



Figure 20. — Change in R t h (normalized) (diamonds), and denormalized 
(squares) with gate width W g . Here n = 10 and 2 l x = 3.26 gm. The 
values from Darwish are 6.3, 9.6, and 1 1 .4 at W g = 40, 200, and 
400 gm, respectively. 


The next comparison considers the work from Nuttinck et al. (ref. 19). This is one of the few papers 
wherein the variation with the number of fingers has been shown. The results are those from a simulation, 
thus we must proceed with some caution. First of all, some of the dimensions, values used for thermal 
conductivity, etc., were not given, but from study of other publications by the authors, we have been able 
to develop a reasonable interpretation of their results. In figure 21 the dimensions of the device are shown 
schematically. The configuration is highly asymmetrical as the source-gate separation is only 0. 1 pm, 
whereas the gate-drain spacing is 3 pm. The thermal conductivity for SiC at 300 K was assumed to be 3.9 
(from one of their other publications). The remaining ones were chosen to be typical values. The 
simulated channel temperature versus the number of fingers is given in figure 22. Notice the temperature 
increases linearly as n varies from 1 to 8. Flowever, the value at n=\6 falls well below an extrapolated 
value from the line developed for small n. Apparently the temperature given is that of the center finger, 
and the authors state that above 10 to 12 fingers, additional fingers do not lead to a significant 
modification of the horizontal heat flow. After 14 fingers, a plateau is reached. Apparently this plateau is 
the value that the one-shot calculation should produce. We choose n = 6 for the one-shot calculation. At 
that point the channel temperature is 435 K and the thermal conductivity of GaN is calculated to be 1.65. 
The thermal resistance for the GaN layer is found to be 41.3 K/W. For the SiC layer we find it to be 
100.6 K/W. The ratio of thermal resistances is then 0.41. From the transient simulation reported in 
(ref. 69), that ratio was 0.6. Using the total thermal resistance of 141.9 K/W, and the dissipated power of 
0.625 W/finger, the calculated channel temperature is (lower value method) 41.3 x 0.625 + 100.6 x 6 x 
0.625 + 20 = 423 °C. Then for «=12 we have 800 °C. These are very high compared with the simulations. 
The GRC-1 method gives thermal resistances of 66.6 for GaN and 34 to 72 for SiC as n varies from 1 to 
16. Figure 23 summarizes these calculations. The prediction of the channel temperature for the GRC-1 
method is also shown as the line in figure 22. Both calculations do not correlate well with the simulated 
curve. Flowever, we can get complete correlation with the GRC- 1 method in the following way. Starting 
with the single finger case we obtain the footprint at the GaN/SiC interface to be 8.35x 127 pm 2 . Then 
with the data set obtained from their publications, we calculate R,i, for the SiC as follows: 

21 x = 8.35 pm 
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2l v = 127 pm 
2 ,L X = 590 pm 
2 L y = 278 pm 
W =330 |im 
p J = 3.9/2.4= 1.625 
Y<? = 15.2 
7, = 0.47. 




Figure 21. — Top and cross-sectional views of device 
simulated by Nuttinck et al. (ref. 19). 
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500 


Simulated and calculated channel temperature 



Number of fingers 


Figure 22. — Channel temperature of simulated device. The lower curve is derived 
from Nuttinck’s simulation. The straight line is the GRC-1 prediction, while the 
single point (n = 6, T= 435) is the one-shot prediction. 


GaN >41.3 


• 66.6 


■ 66.6 


SiC <100.6 <34-72 <26.2 

) n = 1,16 ^ 

(a) ' (b) (c) 

Figure 23. — (a) Is the resistance stack calculation by the 
one-shot procedure, (b) Is that for the GRC-1 method, 
(c) Masana’s method, neglecting cross-talk in SiC. 


We obtain the thermal resistance for SiC to be 26.2 °C/W. The temperature determined by the 
simulator was 78 °C, and subtracting the base plate value of 20 °C, yields a net drop across both the GaN 
and SiC layers to be 58 °C. The drop across the SiC is calculated to be 26.2 x 0.625 = 16.4 °C. Therefore 
the drop across the GaN is 41.6 °C. This gives the thermal resistance across the GaN the value 
66.56 °C/W, with corresponding thermal conductivity 1.8 W/(cm-K). This value is well within the range 
for GaN. Now assuming the change in temperature drop across the GaN is constant as the number of 
fingers increases (since the fingers are sufficiently separated to eliminate any cross-talk in the GaN). We 
obtain the following equation for the channel temperature: 


T ch = R lh (n ) (0.625)« + 20 + 41.6. 


Using the simulated temperatures with n, we obtain R,i,(n) in the SiC. From these we calculate the 
required values of thermal conductivity as a function of temperature. The values needed are shown on the 
solid curve (squares) in figure 24. This curve is reasonably close to the average one shown. The four 
values at 300 K obtained in the literature show the large dispersion one may expect. Comparison with 
figure 22, where constant K SiC was used, shows how critical the correct value for k(7) for any layer must 
be in order to match measured or simulated results. As stated earlier, this data on the thermal 
conductivities is not known precisely, which hampers any analytical scheme. 
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Approximations for SiC thermal conductivity 



— • — 3 . 4 ( 300 / 70' 5 

— ■ — Required by Masama's 
equations 

— *— Required by ATLAS 
simulator 

— x— Published value 
— Published value 
—•—Published value 
—i — Published value 


0 

280 320 360 400 440 480 520 560 


Temperature (K) 


Figure 24. — Composite of thermal conductivity of SiC. The multiple values at 
300 K are 2.8, 3.75, 4.7, and 4.85 are those found in the literature. 


VII. 2-D Silvaco — Atlas Simulations (Basic Information) 

This section encompasses the material parameter sets and computational models used in the 2-D 
simulations. The version of ATLAS used was 5.8.3.R. The flags set were Fermi-Dirac statistics, bandgap 
narrowing, field dependent mobility, SRH generation/recombination, and lattice temperature. Impact 
ionization and tunneling were not allowed, the drift-diffusion option was used, and the baseplate 
temperature was 300 K. No traps were included. At the present time, the information on traps is not 
precise, and a wide range of parameters describing them exists. The physical location, capture cross- 
sections, energy levels, and densities are widely dispersed; and the trap signatures depend on many 
process steps, as well as the growth mechanism, and the substrate. Therefore, we felt including them into 
the model at this time would overly complicate an already very complex model. In most cases, the electric 
field strength was below the critical field for ionization, so we don’t thi nk turning off that option caused 
significant errors. The boundary conditions were the default provided by ATLAS, the Fermi level was not 
pinned, as no surface charges, or traps were included. For any parameter not mentioned in the text, the 
default was used. 

The model was developed with internal self-consistency as a primary goal. We attempted to stay 
within the ranges found in the literature whenever possible. However, the ranges are wide, due to the 
large number of parameters available for the devices. We chose values specifically for comparison with 
the experimental results of Kuball et al. (refs. 15 to 18), wherein both measured I-V characteristics and 
channel temperatures for several drain bias voltages were given. Their work was unique in that both 
accurate self-consistent electrical and thermal information was provided. On the other hand, many of the 
key parameters needed to start a simulation were not given. The chip dimensions, 2DEG sheet density, 
and low field mobilities were given. The mobilities were 1560 and 1900 cnr/V-sec for sapphire and SiC 
respectively. A common electron sheet density of 7.5><10 12 cm 2 was specified. The doping levels were 
not specified, so we adjusted them in the simulation to permit best fitting of both current and temperature 
data. 

The device cross-section is shown in figure 25. Starting at the top, it consists of two Al 2 Ga 8 N layers; 
the upper and lower ones are 230 and 50 A thick respectively. Their doping is 2x 10 18 and 1 x 10 15 donors 
per cm 3 , respectively. The highly doped layer serves as the carrier supply source. The field at the 
AlGaN/GaN interface due to both the spontaneous and piezoelectric polarization is modeled by a positive 
sheet charge ct s = lxlO 13 and 0.7xl0 13 charges per cm 2 for GaN grown on sapphire or SiC respectively. 
These values were chosen by trial and error. Observe the magnitude value for SiC is very close to the 


NASA/TM— 2008-215444 


32 



s 

77771 


G 

IZZL 


D 

V77777— 


M = 600 

230 A Alo.2 Ga 0.8 N 

+ + + + + + + + + 

5o|a 

2 deg channel 

2o!a 

Transition 

io|a 

Smoothing 

2 A 


M = 900 bulk GaN \ 


" ai 0.2 Ga 0.8 N 
,, GaN 


jj =1560 sapphire 

| 1500 

1 A 1440 

t 1380 

1320 

1260 

1200 

1140 

1080 

1020 

960 

2|A 900 


r n =2000 SiC 
1900 
1800 
1700 
1600 

>Transition< 1500 
1400 
1300 
1200 
1100 

J V 1000 

900 


Figure 25. — Schematic (not to scale) of simulated devices. 
Figure 26 gives additional information. 


value for the sheet density given by Kuball. This approximation has been used by several authors in the 
literature (see for example (ref. 31)) with apparent success. In principle, the bound charge due to the 
polarization is a distributed effect; but incorporating this with the available options in the code is not 
straightforward. Due to excellent agreement with published data, we think the single sheet model is 
adequate for our purposes. The strain state of the GaN is different, depending on the substrate, thus the 
different sheet values. When GaN is grown by MOCVD it is in compression/tension depending on the 
substrate being either sapphire/SiC respectively. This difference may cause the polarization effects, and 
traps, to be different. Due in part to these differences, our discussion is separated into two parts; one for 
sapphire and one for SiC. 

The sheet charge values were assumed constant with position along the channel for all gate and drain 
voltages and temperatures. Experimentally, n s the negative charge density existing in the channel, that 
partially shields the fixed positive sheet charge density, decreases with increasing temperature. This may 
or may not reflect changes in the sheet charge value. We have neglected this effect. The simulated 
saturated drain current reduced by ~10 and 1 1 percent as the sheet charge was reduced by factors of 5 and 
50 respectively. The sheet charge and the alignment parameter for conduction band edges provided a 
reasonable band shape near the AlGaN/GaN heterojunction. We assume any excess charge, or interface 
dipoles that may be present in the actual device; are adequately modeled by the alignment parameter and 
sheet charge. The position of this charge sheet is important; and we get best results when it is directly at 
the physical heterojunction. In actual devices, the surface passivation is a critical process; as much of the 
dispersive behavior is dependent on the charge state in and around the passivation layer. We have omitted 
any charge build-up on the surface, or in the AlGaN. We have also neglected possible conduction paths 
near the GaN/AIN interface. However, we have assumed the AIN to have a doping level the same as the 
bulk GaN layer. The underlying sapphire or SiC substrates have been assumed completely insulating. 

The 2DEG channel (defined by high mobility and specific variation with temperature) was modeled 
as a layer 20 A thick. Below the channel is a 10 A transition layer. This layer is further partitioned into 
ten 1 A layers; wherein the mobilities transition from the channel to the bulk values. Finally, a 2 A 
smoothing layer with the bulk mobility is below the transition layer. Originally this layer was just a 
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placeholder. The mobility values used in the transition layer are listed in the expanded portion of the 
figure. Below this layer is the large bulk GaN layer. Table 9 provides the dimensions for the source, gate, 
and drain pads and their separations. 


TABLE 9.— DIMENSIONS OF SOURCE. GATE, 
AND DRAIN PADS FOR THE BASIC DEVICE 


Element dimensions, 
|im 

Sapphire 

substrate 

SiC 

substrate 

Source pad 

50 

50 

Gate length 

4 

1 

Drain pad 

50 

50 

S-G spacing 

3 

1 

G-D spacing 

3 

2 

Overhang (semiconductor) 
beyond drain edge 

10 

10 


The simulated device of figure 25 is 1 pm deep (into the page), and the results are scaled to the full 
width of 200 pm (done internally in ATLAS). The minimum grid spacing in the x and v-directions was 
0.1 pm and 1.0 A, respectively. The fine mesh in the v-direction was in and about the channel layer. The 
number of grid points for the sapphire and SiC substrate cases were 19,968 and 15,808. The maximum 
number allowed was 20,000. For our simulations the density of points was the same; the larger number in 
the sapphire case was due to the device just being larger. Some parameter values used are listed in table 10. 


TABLE 10.— PARAMETERS FOR THE PRIMARY LAYERS OF THE HEMTS 


Parameter 

Al 2 Ga.gN 

GaN 

AIN 

SiC 

L„eV 

3.96 

3.4 

6.14 

3.34 

N c , cm -3 

2.07xl0 ls 

1.07xl0 18 



N v , cm -3 

1.16xl0 19 

1.16xl0 19 



E d , eV 

0.044 

0.044 


0.0518 

E a , eV 

0.045 

0.045 


0.195 

£,- 

9.5 

9.5 

10.4 

9.66 

N d , cm"' 1 

2xl0 1! ’ upper 
lxlO 15 lower 

l.OxlO 17 

l.OxlO 17 

0.0 

Fsat, cm/s 

1.12xl0 7 

2xl0 7 




The Schottky barrier height was 4.5 eV for all cases, on both substrates. We swept its value from 4 to 
5, and settled on 4.5. The corresponding surface densities n s were 7.4, 7.9, and 8 . 1 3 x ] 0 1 2 cm 2 
respectively. This was performed at the gate center. 

As noted in previous sections, the thermal conductivity values are critical, as the developed 
temperature due to self-heating is very sensitive to the value for k (T). Table 1 1 gives the information 
concerning the k(T) specifications. The thermal conductivities for air above the device, the contact pads, 
and the AlGaN were 2.7- 10A, 2.4, and 0. 13 W/(cm-K), respectively. For the materials with temperature 
dependent k, we used the format, k = l/(aT 2 + bT + c ), and chose 3 points to determine the values for a, 
b, and c. For sapphire the points (k ,T) chosen were (0.486, 363 K), (0.416, 403 K), and (0.40, 445 K). For 
SiC, (3.6, 350 K), (2.6, 400 K), and (1.95, 450 K). Notice these points are plotted in figure 24, and show 
good correspondence with the average curve developed from data in the literature. For the GaN layer, the 
points were derived from the analytical formula 1.6 (300/T) 1 ' 4 . In all cases the sapphire and SiC were 
assumed undoped. In one run we assumed the SiC was doped at 10 15 cm 3 , and the drain current at Vd = 
20 V increased by 3.75 and 5 mA, over baseline values of 27.5 and 57 mA, respectively. These were for 
the gate at -4 and -2 V, respectively. 
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The specific contact resistivities to the GaN layer on the sapphire and SiC substrate cases were 
5.0* 10" 4 and 1.0-10 4 Q-cm 2 , respectively. Variation of these numbers caused a strong effect. The values 
seem to control the linear slope for low drain voltages to a large degree. 


TABLE 1 1.— THERMAL CONDUCTIVITIES (W/CM-K) FOR THE PRIMARY LAYERS IN THE STRUCTURE 



AlGaN 

GaN 

AIN 

Sapphire 
(see text) 

SiC 

(see text) 

K 

0.13 

1.6(300/7)' 4 

0.1 

U(aT 1 + bT+c) 

U(aT 1 + bT+c) 


The carrier mobilities were chosen as follows. In both the AlGaN and GaN layers the dependence on 
the electric field was 


p(e)=\ 


1 + 


Vo E ' 2 

v u sat 


(49) 


In the AlGaN the coefficient p 0 has the form 


p o (r)=p o (300K) 


" 300 " 

V T , 


TMUN 


(50) 


where the value at 300 K was assumed to be 600 cm 2 A" -sec, and TMUN = 1 .5. In all subsequent 
discussion, the units of the mobility will be omitted for brevity. In the GaN channel and transition layers, 
TMUN was 3.1 and 1.0, respectively. These values were assumed for both substrate cases. The 
dependence on doping and temperature in the GaN bulk was modeled using the Albrecht et al. (ref. 70) 
approach. The equation is 


1 

W) 


100 v K a 300 


,-1.5 


,1.5 


+ b\ 


, 300, 


exp 


1065 


-1 


(51) 


= 64.6 


300 


The value for p„( 7) is a fitting parameter, and p„ (300 K) = 900 gives an excellent fit for sapphire. 
The fitting coefficients ( a , b, and c in the above equation) were determined for our assumed doping of 
10 15 /cm 3 . The values used were a = -2.71T0 -2 , b = 2.64-10” 3 , c = — 1 .34- 10 2 . The room temperature bulk 
and 2DEG channel mobility values on sapphire and SiC substrates were 900, 1560 and 900, 2000 cm 2 /V 
respectively. These values correspond to many theoretical predictions and several measured values 
(refs. 71 and 72). Figure 26 takes a more expanded view of the device, and in particular showing the 
contact depths from the heterojunction plane. 

The band diagrams for both substrates are depicted in figures 27 and 28. 


NASA/TM— 2008-215444 


35 




Band diagram — Sapphire 



Figure 27. — Energy band diagram for sapphire substrate case. 
Band diagram — SiC 
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Figure 28. — Energy band diagram for the SiC substrate case. 

In theory the conduction band energy is composed of many parts. For our purposes in this section, we 
will decompose the total band edge energy into two parts; one due to the electric potential, and the second 
due to “other” terms. The “other” terms will be treated in detail in the Discussion section. By definition, 
the electric potential is a continuous function at all points; including interfaces. Thus the potential energy 
associated with the electric potential is continuous everywhere. At a heterojunction, a large change in the 
“other” part occurs over a very small region (say 5 A), which is modeled as a jump in the band edge. 
ATLAS uses the following method to incorporate the energy jump. The magnitude of the jump (the step 
height) is the product of the difference in band gaps of the AlGaN and GaN, and the align parameter. We 
chose the parameter to be 0.8, which gives a step of 0.45 eV. We specified the electron affinity of the 
AlGaN to be 3.82 eV, but after adjustment in the program, its effective value was changed to 2.83 eV. 

The adjusted electron affinity for the GaN was 3.28 eV. The energy band diagram is then constructed as 
follows. The electrostatic potential is determined throughout the structure; then the step of 0.45 V is 
added at the plane of the heterojunction. This step exists over 3 mesh points; one point is at the edge of 
the AlGaN, one in the interface, and one at the edge of the GaN. The function is then multiplied by -1 and 
given the units of electron-volts. Notice the “other” parts are incoiporated solely at the interface, and the 
remainder follows the shape (negated) of the electric potential. 

The sheet charge is placed at the interface, and perturbs the potential in this region. This causes the 
resulting notch depth to be near 0.38 and 0.28 eV below the Fermi level for the sapphire and SiC cases, 
respectively. The size of the step is held constant by the program. From experiments, the measured step 
for our case should be near 0.45 eV. No fixed charge is associated (in Poisson’s equation) with the step, 
as its origin is due to “other” parts. Recall the electric potential in Poisson’s equation is associated with 
both free and bound charge densities, and the charges that produce the step are not included. The effect of 
these unspecified charges is incorporated in the align parameter. Actually the step is due to the differences 
in the bandgaps between AlGaN and GaN. The bandgaps are due to the corresponding lattice constants 
and periodic crystal potentials in the two materials. The sheet charge causes a step in the normal 
component of the electric displacement as expected. The sheet charge is not visible in plots due to either it 
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not being included as an option to be plotted, (since it only exists along a single mesh line) or hidden by 
the line defining the interface. However, its effects are seen in the potential shape and all other applicable 
parameters. When the charge sheet was placed at various planes in the AlGaN, small notches appeared 
about them; which we considered to be unsatisfactory. This is apparently due to the software initially 
placing neutralizing free charge at adjacent mesh points to keep the entire system charge neutral. When 
this charge is placed in the notch, most of it apparently stays there or in nearby regions. When the sheet 
was moved to different locations, most of the free charge moved away from the plane, and left some 
residual near the plane. That residual along with the sheet charge caused a small well about the plane of 
the sheet. This is also problematic when a negative sheet charge is placed on the top surface of the 
AlGaN; holes are immediately inserted and must be dealt with. There is no easy way of handling the 
holes since they should not be present in the material. The bandgap is too large for thermal generation. So 
once placed in the system, they cannot combine, so what should one do to address their artificial 
presence? One might place some (always ionized) compensating acceptors on the surface to hold them. 
However, many numerical simulations we reviewed in the literature did not place any charge on the top 
surface; and apparently achieved good results. Therefore, we chose not to use more than one charge sheet. 


VIII. Sapphire Substrate Results 

The simulated (solid) and measured (dotted) I-V characteristics for a HEMT on a sapphire substrate 
are shown in figure 29. These data were measured by Kuball et al. (ref. 15), and we assume no trapping 
effects were present. Thus the I-V should be reasonably independent of sweep rate, direction of gate step 
voltage, and we assume the curves were developed in the thermal steady-state. The droop is assumed to 
be solely due to self-heating. In most practical devices the droop can be due to both self-heating and 
trapped charge in various regions of the device. We have neglected charge trapping, and apparently these 
devices have minimal trapping effects. When the lattice temperature flag was turned off, the curves did 
not droop. 

The V g = 0 curve is the primary one we will focus upon, as the measured temperatures were 
determined for this gate voltage. In general, the match between simulation and measurement is excellent. 
When approaching pinch-off, we see the conduction in the bulk is slightly too high. There are many 


AIGaN/GaN/Sapphire HEMT 



Figure 29. — Simulated (solid) and measured (dotted) I-V 
characteristics for an AIGaN/GaN HEMT on a sapphire 
substrate for gate voltages of 0, -2, and -4 V. 
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Lattice heating effect 
Sapphire 



Figure 30. — The change in drain current when lattice heating 
was turned off or on (for the V g = 0 case only). 

changes we can put forth to remedy this; we could reduce the mobility and/or saturation velocity and/or 
the doping in the lower region of the GaN where the major portion of the current exists near pinch-off. 
Figure 30 shows the changes in the I-V curves (only the effect on the V g = 0 curve) for lattice heating 
turned off or on. 

When the lattice heating was turned off and zero contact resistance was also in effect, the V g = 0 curve 
reaches slightly above 100 mA at 20 V. Figure 31 gives the comparison between the simulation and 
measured temperatures. The comparisons are made at drain voltages of 4, 8, 12, and 20 V. The vertical lines 
with dots are just to aid the eye. The simulated temperatures at the gate edge are below the curve, and the 
ranges of temperatures from Kuball’s paper are given above. The simulated values tend to the low sides of 
the ranges at low drain voltages, and then to the high side near 20 V. Only at 12 V is the simulation 
comfortably in the range. The measured data ranges are the same as those given in table 2 in section V. 
Observe the ranges are not precise, and only a visual judgment of the measured points was used. 

Figures 32 and 33 give the temperature variation (drain at 20 V, gate at 0 V) along a cutline through 
the middle of the gate to the bottom heat sink. 

The contour plot in figure 32 shows that the isotherms become horizontal near a depth of 50 pm. The 
peak temperature is 197 °C, and that at 50 pm is 137 °C. The right panel shows the nearly linear drop in 
temperature through the sapphire. Figure 33 gives the detail near the surface. The drops across the GaN 
and AIN regions are about 2.1° and 0.3° respectively. Figure 34 shows the drop across the device for 
increasing drain voltage (unless stated otherwise, the gate was always at zero). We can generalize; a 
nearly exponential drop occurs for the first 50 pm, with the remainder being linear. The linear drop occurs 
when the isotherms are horizontal (as expected). Figure 35 shows the peaking behavior at the gate edge 
(on the drain side) as one moves horizontally from source to drain. The spike occurs at the gate edge in all 
cases. Kuball’s laser spot size was about 1 pm 2 , and his longitudinal steps were between 0.2 and 2 pm, so 
his measurements should capture the correct average temperature. 
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Lattice temperature (°C) 
Sapphire 



Drain voltage (V) 

Figure 31. — The measured range of channel temperature (above 
curve) and the simulated values (below) at drain voltages of 4, 
8, 12, and 20 V. 


Device temperature 
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Figure 32. — The temperature contour plot and value along a vertical cutline through the center of the gate for the 
drain at 20 V. 
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Detail of temperature profile 
Sapphire (V d = 20 V) 




Figure 33. — The detail of the temperature in the AIGaN/GaN and AIN regions (left panel) along with the variation 
through the sapphire (right panel). 


Figure 36 gives temperature changes as one moves from the gate towards the drain. The three curves 
are values at the channel, middle and bottom of the GaN layer, respectively. At the gate edge the values 
are 482, 478.5, and 475 K; which is only a difference of 7°. Therefore averaging over the depth should 
not produce significant error. 
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Temperature profiles through device 
Sapphire 



Figure 34. — Temperature drops across the sapphire substrate for 
drain-source voltages of 4, 12, and 20 V. 
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Temperature profiles in GaN 
Sapphire 



Figure 35. — Variation of temperature as one moves horizontally 
through the center of the GaN layer. 
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Temperature profiles between gate and drain 
Sapphire 



Figure 36. — Temperatures between the gate and drain edges at the top, 
middle, and bottom of the GaN layer. 
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Temperature profile between source and drain 
Sapphire 




Figure 37. — A thermal contour plot and the temperature from source to drain at the middle of the GaN layer. The 1 A 
thick transition layers are shown in the contour plot. 


The horizontal cutline in figure 37 is at the center of the channel in the GaN layer, and it reveals the 
change in temperature from the source to drain contacts. The small kink in the curve in the panel on the 
right is at the gate edge near the source. The source and drain edges are at 179 and 187 °C, and the edges 
of the gate are at 190 and 209 °C. 
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Table 12 shows the drain current at Vd = 20 V, at zero gate bias, for two different variations on 
temperature of the thermal conductivity. Observe the current differs by 3 mA while the temperatures 
differ by 32°. Smaller differences are seen for the gate at -2 V. 


TABLE 12.— CHANGE IN PEAK TEMPERATURE AND DRAIN 
CURRENT FOR TWO DIFFERENT APPROXIMATIONS FOR 
THE THERMAL CONDUCTIVITY FOR SAPPHIRE 


K (T) 

Id( 20 V), 

T, 


mA 

K 

1 IT 

31 

482 

l/(af- + bT+c) 

34 

450 


Figure 38 shows the temperatures at the gate edge (at the surface) and at vertical planes in the gate- 
drain access region of the GaN layer versus the drain voltage. 

Figure 39 reveals the source and gate edges track each other with drain voltage. Figure 40 is a 
potential contour plot; the steps are in 3 V. 


Peak temperatures at selected depths 



Figure 38. — Gate edge temperatures and that in the gate 
drain access regions as the drain voltage is increased. 


NASA/TM- 


-2008-215444 


47 


250 


Temperature in GaN 



Figure 39. — The temperatures at the gate and source edges with increasing 
drain voltage. 


Potential plot at Vd = 20 V 
Sapphire 
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Figure 40. — Constant potential contours (3 V steps), 
showing the source, gate, and drain pads from left 
to right. The gate is at 0 V. 
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Potential (V) 


Figure 41 shows the potential at a vertical outline through the gate on the left, and along a horizontal 
outline at a depth of 0.0501 pm. 


Potential under gate and in channel 
Sapphire (V d = 20 V) 




Figure 41 . — Variation of potential along a vertical outline near the gate center (left panel), and along a horizontal 
outline along the channel at a depth of 0.0501 pm for drain (right panel) at 20V. 
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Electric field in AIGaN and channel 
Sapphire (V d = 20 V) 




Figure 42. — The electric field in the AIGaN from source to drain (left panel). That in the GaN is shown on the right; 
Observe the field drops two orders of magnitude at the gate edge. 


The electric field magnitude from the source to drain in the AIGaN and GaN layers is depicted in 
figure 42. The left panel shows the peaking at the gate edge in the AIGaN buffer region. Observe it 
approaches 3 MV/cm, which is close to the avalanche value. The literature states the avalanche field is 
somewhere between 3 and 6 MV/cm. The right panel shows the much lower value in the channel. 

Observe the notch at the gate edge where the field drops to 1 .8x 1 0 4 V/cm. Figure 43 shows the magnitude 
of the electric field near the channel, from source to drain for drain voltages of 20 (left) and 8 (right) V. 
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Electric field near channel 
Sapphire 




Mesh coordinate Mesh coordinate 

Figure 43. — The variation of the magnitude of the electric field near the channel for 20 V (left panel) and 8 V (right 
panel) on the drain. 


The magnitude of the velocity field at two horizontal outlines is illustrated in figure 44. The left panel 
is at a plane below the channel, while that on the right is at a deeper plane. 

Figures 45 through 48 demonstrate the carrier concentration variation along specific outlines. 
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Electron velocity 


Velocity field below channel 
Sapphire 
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Cutline below channel 
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Cutline further below channel 


Figure 44. — The velocity field (cm/s) at two depths below the channel. 
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Carrier densities under gate edge 
Sapphire 
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Distance from interface 

Figure 45. — The distribution of carriers along a vertical cutline at 
the gate edge through GaN at Vd = 12 V. 


NASA/TM— 2008-215444 


53 






- 0.2 


Carrier density in GaN and AIGaN 
Sapphire 



pm 



Figure 47. — Contour plot of the carrier density. The right panel shows the concentration along a horizontal outline 
through the AIGaN at a depth of 0.030 pm. Notice the region under the gate is completely depleted. 
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Current density in AIGaN 
Sapphire 




pm pm 

Figure 49. — The current density in the AIGaN layer versus horizontal position. Cutline at 0.035 pm. 


Figure 50 demonstrates the interesting behavior of current density variation with depth into the GaN. 
The current drops in the lower half, but then picks back up near the AIN interface. Notice the absence of 
current at the heterojunction, which means a portion of the “channel” is completely depleted. Therefore in 
the saturation region, the current is supported by conduction solely in the GaN layer over horizontal 
intervals encompassing the gate edge plane. Further down the channel (say midway to the drain) the 
current in the channel picks back up. A spike of current then exists at the heterojunction. 
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Current density under gate edge 
Sapphire 



pm 



Figure 50. — Contour plot of the current density. The right panel shows the current density along the vertical outline at 
the plane of the gate edge 


Figure 5 1 has an interesting contour plot of the current density. The light regions indicate larger 
values. The plot on the right is along the cutline shown just below the channel. Notice the deep notch at 
the plane of the gate edge. 


NASA/TM— 2008-215444 


58 



Current density under along channel 
Sapphire 



pm 



Figure 51 . — Current density along a horizontal outline slightly below the channel. The notch occurs at the gate edge. 


Figure 52 demonstrates the heat generation profiles at horizontal cutlines. 
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Sapphire 




Figure 52. — The heat generation profile at two different horizontal outlines. The drain is at 20 V. 


Figure 53 is a close-up view of the heat generation rate along a horizontal outline in the channel area. 
Observe the heat generation at the source and drain contacts. The heat generation is composed of several 
terms; the dominant one being the dot product between J and E, the other components are assumed small. 
The notches occur in the portions where the current is very low. 
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Heat generation rate at V d = 20 V 
Sapphire 
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Figure 53. — Contour panel and outline along the channel of the heat generation rate (W/cm 3 ) for 20 V on the drain. 
Notice the heat generated at the ohmic contacts. 


Figures 54 through 57 illustrate the heat generation rate as the voltage is increased from 4, 8, 12, and 
20 V. The shapes of the plots leads one to speculate that a combination of point sources (for the spikes) 
along with uniform regions (averages of the remainder), may yield better results in analytical models. 
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Heat generation rate (W/cm3) 





Heat generation (W/cm3) 







In general there are four spikes in g(x,y), the heat generation rate. They are at the contacts and edges 
of the gate. It is interesting that g(x,y ) is near zero at the singularity plane due to the current (the current 
density is as low as 3 A/cm 2 ) being near zero there. In spite of this, however; the temperature peaks in this 
region. We close this section by stating that the channel current is disrupted in an interval about the gate 
edge. For outlines not exactly at the heterojunction, the current is nearly zero. Much of the current is 
carried in the upper portion of the bulk GaN layer in the saturation region. Thus the detailed properties of 
the channel are not especially important for the device current. The properties of the upper portion of the 
GaN layer are important for current. It is also interesting that a portion of the current exists near the GaN- 
A1N plane. This appears to be due to the geometry of the device; that is, a favorable horizontal component 
of the electric field. In some cases the channel is re-established near the drain pad. If both the gate-drain 
spacing is small, and the voltage difference is large, then the channel may not re-establish itself. 

IX. SiC Results 

The measured and simulated I-V characteristics of the SiC substrate device are given in figure 58. 

The agreement is not as good as that for the sapphire case. The calculated pinch-off currents don’t reach 
the measured values for low drain voltages and the curves keep on rising at the higher voltages. There are 
several reasons that explain this discrepancy. First of all, the initial SiC device simulation used all of the 
parameters found earlier for the sapphire case; the only changes incoiporated were the source-gate-drain 
dimensions and the thermal conductivity value. The results were not good. We varied many of the same 
parameters that we found useful in the sapphire case, but could not get all of the three gate voltage curves 
to simultaneously match the experimental ones, and agree with temperatures. Only after altering the 
supply layer density in the AlGaN (something we did not do for the sapphire), did we start to get much 
better correlation. Unfortunately, time ran out for this project; so figure 58 is given as our best effort. The 
increase in doping in the supply layer by a factor of two brings the V g = -2 curve up to the measured 
value. This fact shows the sensitivity of the parameter set needed for good agreement throughout the I-V 
plane. In particular, the supply layer is critical because electrons cannot come from anyplace else to form 
the conduction path. Unlike some experimental devices wherein just the unintentional residual doping 
supplies enough carriers to establish drain current; unless the supply layer density is sufficient, the 
simulation currents will all be low. Figure 59 shows the sensitivity of the drain current as the supply layer 
doping is increased by a factor of two. The gate equal to zero and -4 curves were determined for doping 
of 2x 10 18 cm' 3 , while the three central curves were for the gate at -2 and the changing doping values. 

This increase shows the short- fall in the previous plot was due to the lack of carriers available in the low 
drain voltage regime. In other words, the increased doping prevents the gate potential from pushing too 
many carriers out of the channel region. We observed, but do not show, that the saturated drain current for 
two values of channel mobility are rather close, which shows the I-V is more sensitive to the supply 
doping value than the 2DEG mobility. The maximum temperature in the sapphire case was about 1 80 °C, 
and that for the SiC was closer to 120 °C. The current in the SiC case was about twice that of the sapphire 
case. Thus the SiC case was both cooler and with stronger fields. These factors demonstrate the 
challenges presented in higher power devices, smaller geometry devices; and the corresponding need for 
closer scrutiny. We choose, however, to ignore these factors at the present time, and leave them for future 
investigations. All of these factors must be considered; and the nonlinearities are unfortunately not known 
precisely, which makes the current very sensitive to parameters. Thus one must have “everything right” 
for the simulations to agree with experiment. Not to mention the fact that traps may be present, which we 
have neglected from the outset. Even with these caveats, we think the agreement with measurements is 
very satisfactory. 
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— * — Vg = 0 (simulation) 
—a — Vg = -2 (simulation) 



Figure 58. — Static l-V characteristics for SiC substrate 
HEMT. The legend indicates the simulation versus 
experimental curves. 


SiC l-V Curves 



Figure 59. — The increase in the drain current for V g =-2 as 
the doping density of the supply layer is increased by a 
factor of two. 
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Figure 60 gives an expanded view of the temperatures near 20 V. The difference in temperature 
between the gate edge and the bottom of the GaN layer for this case is about 15°. The results of (ref. 68) 
show consistently higher temperatures when just the surface temperature is measured. Thus we have some 
leeway when matching-up simulation and measurement. The average is not well defined, so we need only 
to be within say 20° to declare a good match. In other words, is the measured temperature the average of 
the GaN layer, or the peak value on the top surface of the AlGaN? A final remark is that the temperatures 
in horizontal planes in the GaN show the gate edge and source edge track each other. At low drain values 
they are nearly the same, and at 20 V the source is about 18° below the gate. 

The electron concentration along a vertical cutline through the center of the gate is given in figure 6 1 . 
Observe it is climbing from essentially 0 to 10 14 cm at the heterojunction. It peaks near 8 x 1 0 1 x cm 3 in 
the notch. Figure 62 depicts the temperature drop across the AlGaN, GaN, and AIN layers. In most cases 
the drop across the GaN layer is around 22°, while that across the AIN is near 3°. The AlGaN is only a 
fraction of a degree hotter than the GaN. Figure 63 gives the thermal profile through the entire device. 
Similar to the sapphire case, an almost linear drop occurs from around 50 pm depth down to the heat si nk . 
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Figure 60. — Temperatures at the gate edge and in the gate-drain 
access region as the drain voltage increases. 
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Temperature variation in AIGaN, GaN, AIN 
SiC (V d =20 V) 



Distance from interface 

Figure 62. — Temperature profile in the AIGaN, GaN, and AIN 
regions. The drop across the GaN is typically in the 21° 
to 24° range for most conditions, such as mobility 
variations and doping levels. 


The two line plots in figure 67 are horizontal cuts slightly below the channel for 4 V on the drain. 

Both show heat generation, and that on the right is at a slightly lower depth. Observe the two spikes that 
exist near the gate edges in the access regions. The drain access region is hotter by nearly an order of 
magnitude. Figure 68 gives vertical cuts at x=52 and 54 pm, respectively. On the left the cutline is defined 
by 0.476<v<0.094, while on the right it is 0.046<v<0. 131. Atx = 52 (gate edge), g(x,y) extends 0.01 pm 
into the material. At x - 54 (drain edge) we have the ohmic spike, a drop in the AIGaN region, then the 
spike at the heterojunction. 
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Device temperature 
SiC (V d =20 V) 



Figure 63. — The temperature profile through the entire device; from 
gate to heatsink. The drop in the SiC is essentially linear for the 
lower 300 pm of the substrate. 


Thermal contours at 4 V on the drain (V g = 0 V) are illustrated in figure 64. The steps in shading are 
only a few degrees apart as most of the device is near room temperature. Table 13 demonstrates the 
temperature variation at arbitrary points along the upper portion of the device. 


TABLE 13.— SIMULATED TEMPERATURES AT VARIOUS 
POINTS IN THE UPPER REGION OF THE HEMT 


X 

y 

Temperature, 

K 

0.84 

3.6 

307 

29.9 

3.64 

307 

51 

2.8 

309 

76 

2.8 

307 

100.5 

3.6 

307 

113.4 

2.8 

307 
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The coordinates correspond to the contour plot in figure 64. Notice the GaN layer extends down to 
1.2 pm, with the remainder below being SiC. All of the sampled points were chosen in the SiC. The first 
and last entries (not shown in fig. 64) are at the extreme left and right sides of the unit. The baseplate is 
kept at 300 K. 

The heat generation contour plots are exhibited in figures 65 and 66. The drain is at 4 V, and notice 
the heat generated at the ohmic contacts is comparable to that generated in the channel. The small bump 
at the AIN interface is due to a non-negligible amount of current existing there. A closer view of the 
channel region is shown in figure 66; the behavior at the source pad is interesting. 


Temperature contours 
SiC (V d = 4 V) 
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Figure 64. — Constant temperature contours centered on the gate 
edge for bias ( V g = 0 V, V d = 4 V). Observe the interesting 
focusing effect of the low conductivity AIN layer. 
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Heat generation 
SiC (V d = 4 V) 
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Figure 65. — Heat generation at 4 V on the drain. The heat 

generated at the contacts is comparable to that developed in 
the channel. 
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Detail of heat generation 
SiC (t/d = 4 V) 



Figure 66. — Tighter view of the heat generation function 
g{x,y). The behavior near the source pad is interesting; 
the temperature distribution of the air between the 
source and gate metal pads is similar to that in the 
solid AIGaN. 
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Heat generation under gate and drain edge 
SiC (\4/ = 4 V) 




Distance from interface Distance from interface 

Figure 68. — Vertical cuts showing g(x,y ) versus y at planes corresponding to the gate and drain edges (left and right 
panels, respectively). On the drain edge the first spike is due to the ohmic contact. 


Figure 69 demonstrates the mobility and electric field magnitude for 5 V on the drain. The horizontal 
cutline is at y = 0.076 pm. The field peaks at the source, gate, and drain edges respectively. Figure 70 
shows the entire simulation area and elucidates the temperatures at selected points in the structure for the 
bias point (V g = 0V, V d = 5 V). The chip edges on the top surface are at 3 10 K and the base is at 300 K. 
This is in contrast to Kuball’s assumption that the baseplate rises to nearly the value found at the chip 
edges. 
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Mobility and electric field in channel 
SiC (\4/ = 5 V) 




Figure 69. — The mobility and absolute value of the electric field along a horizontal outline near the channel. 


Figure 71 is a contour plot (Vj = 5 V) with nearly the same scale for both axes. Therefore, the 
isotherms have nearly their true shape. The drop is only about 4° from the point in the light region 
about the gate and the lowest contour. Again the isotherms are nearly horizontal at the depth of 50 pm. 
Figure 72 is a similar plot for the drain increased to 8 V. 

Figure 73 shows the magnitude of the current density in the AlGaN/GaN regions. Table 14 gives the 
values at the specified points. 
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Temperature contours through entire device 
SiC (V d = 5 V) 



Figure 70. — View of the entire device showing the 

temperatures at specific points. The gate is at 0 and 
the drain is at 5 V. The drain current is about 98 mA. 


Temperature contours 
SiC (\4/= 5 V) 



Figure 71 . — Contour plot at 5 V with nearly identical 
scales for vertical and horizontal axes. The true 
shapes of the isotherms are seen clearly. Observe 
that the temperature drops about 4° from top to 
bottom. 
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Temperature contours 
SiC(\4,= 8V) 



Figure 72. — True shape of contours for 8 V on the drain. 
Here, about 15° is dropped from channel to bottom 
(about 106 pm away). 


Total current density 
SiC (V d = 8 V) 
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Figure 73. — Current density contours for Vd = 8 V. The 
geometry of the device causes the small portion at the 
GaN/AIN interface. 
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TABLE 14.— THE MAGNITUDE OF THE CURRENT DENSITY IN 
THE ACTIVE DEVICE REGION FOR 8 V ON THE DRAIN 


Point 

X, 

jim 

y> 

|xm 

\J(x-y)\, 

A/cm 2 

A 

52.0 

0.122 

148,702 

B 

52.0 

0.283 

23,624 

C 

52.0 

0.979 

2,774 

D 

50.4 

0.034 

5,222,558 

E 

53.2 

0.052 

7,949,000 


Figure 74 is the x-component of the current density, for the same bias point (V g = 0V, V d =8 V). 
Table 15 gives the values at points A, B, and C. 

Figure 75 accompanies the two previous figures, as it displays the y-componcnt of the current density. 
Figure 76 is a horizontal cutline plot showing J y for V d = 8 V. The two spikes occur at the gate edge 
and the drain contact. The small dip at 1 pm is on the source-side gate edge. 


TABLE 15.— THE x-COMPONENT OF THE CURRENT DENSITY 
TIN THE ACTIVE REGION AT 8 V ON THE DRAIN 


Point 

X, 

(dm 

y, 

\im 

J(x,y), 

A/cm 2 

A 

52.0 

0.204 

50,608 

B 

52.0 

0.531 

6,279 

C 

51.7 

1.23 

8,079 


J x current density 
SiC (Vd = 8 V) 

■ 0.2 
0 

0.2 
0.4 
0.6 
0.8 
1 

1.2 
1.4 


49 50 51 52 53 54 55 

|xm 


Figure 74. — Contour plot of J x at Vd = 8 V. 
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Figure 75. — Contour plot of J y at V d = 8 V. Notice the absence (or change in opposite 
direction) at the gate edge. The current is pushed down into the GaN layer at the 
gate edge. 
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Jy current density below the channel 
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Figure 76. — Horizontal outline demonstrating the behavior of 
Jy slightly below the channel. 


The contour panel in figure 77 is that of J x , while the outline plot on the right encompasses the 
heterojunction (see outline 1 in the left panel). Table 16 gives the values at the indicated points. Point A in 
the AlGaN has considerable current, while point B (in the undoped region of the AlGaN) has almost 
none. Point C is at the heterojunction and the current density peaks there. Point D is just below the 
transition layers, and is the start of the bulk GaN region. Finally point E is deep into the GaN layer. 


TABLE 16— THE x-COMPONENT OF THE CURRENT 
DENSITY J FOR THE FIVE POINTS IN FIGURE 77 
UNDER THE CENTER OF THE GATE (V d = 8 V). 


Point J x Current Density (A/cm 2 ) 

A 2.6x1 (f 

B -238 

C 1.75x10" 

D -1.24x10 s 

E -2106 


The conduction band edge at three vertical planes is illustrated in figure 78. The drain is at 8 V. The 
coordinates of the planes and ends of the cutlines are: x = 50.5, -0.009 < y < -0.095 (midway between 
source and gate); x = 51.5, -0.024 <y< -0.1 (mid-gate); x = 53, -0.009 <y< -0.093 (midway between 
gate and drain). The notch remains all along the channel for this bias point (V g = 0 V, = 8 V). 
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J x current density between gate and drain 
SiC(t/ d =8V) 




Figure 77. — Conduction current contours ( J x ), and its values along cutline 1 . Observe the current is contained in the 
channel region (the plane of the heterojunction and about 0.003 pm (30 A). 
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Conduction Band Profiles 
SiC (t/d = 8 V) 
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Figure 78. — Conduction band edge profiles at mid source-gate access, center of gate, and mid gate-drain access 
region. 

The heat generation plot for 8 V on the drain is shown in figure 79. Notice values in contours that 
contain the points A and B differ by four orders of magnitude. This demonstrates the large heat spike at 
the gate edge, and why an equivalent heat source close to the S-D distance causes the analytical 
expression of Masana to agree with simulations and measurements. 
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Figures 80 and 81 are horizontal cutline plots of g(x,y ) at y = 0.0474 for 8 and 20 V on the drain. The 
ranges of x are: 50 < x < 54.2, and 49.9 < x < 54.2 respectively. Notice the spike (at 8 V) at the source 
contact is the largest at this depth with the gate — drain access being the next largest. Flowever, at 20 V the 
drain spike starts to develop, but is still an order of magnitude below the source’s value. The values in the 
access regions around the gate are nearly the same in both cases; about 4x 10 9 W/cm 3 . 

Figures 82 to 85 are thermal plots for Vd = 20 V. The axes are nearly the same in figure 82 to show 
the true shapes of the isotherms. A closer view is presented in figure 83 to demonstrate the high 
temperature column beneath the gate edge. Figure 84 displays the entire device; notice the chip edges are 
about 58° above the heatsink. The heatsink was held at 300 K, and its top surface was determined to be 
304 K. Apparently the 25 pm thick copper bottom surface is held at 300 K, and the boundary condition 
there is considered to be Dirichlet. The area of the copper surface is 113x200 pm 2 and its thic kn ess is 
25 pm. Using the general expression for thermal resistance R,h = //(kA), and using k = 4 °C/W we obtain 
R,h = 2.765 °C/W. For our simulation we obtained about 91 mA at 20 V, so the dissipation was 1.82 W. 
This yields a temperature drop of 5 °C; which agrees within 1° of the simulation. A cutline through the 
gate edge down into the SiC is depicted in figure 85. The channel is at 430 K and the GaN/AIN interface 
is at 400 K. The drop across the AIN is about 3°. 


Heat generation rate 
SiC (Vd = 8 V) 



Figure 79. — Heat generation around the channel and into the 
GaN layer. 
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Heat generation rate in channel at V d = 8 V 
SiC 
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Figure 80. — Heat generation along horizontal outline in the 
channel area at 8 V. 
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Heat generation rate in channel at = 20 V 
SiC 


8e+09 

4e+09 


-4e+09 

-8e+09 

-1.2e+10 

-1.6e+10 


Figure 81. — Heat generation along same outline as figure 80, 
but now at 20 V. 
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Figure 82. — True shape of the isotherms for 20 V. 
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Temperature contours under gate 
SiC (Vd = 20 V) 
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Figure 83. — Contours of temperature centered on the 
plane of the gate edge for 20 V bias, the gate is at 
0 V. 


Temperature contours for entire device 
SiC (Vd = 20 V, /d = 91 mA) 



pm 


Figure 84. — Contour plot of the entire device. The 
temperatures at the upper edges are around 60° 
hotter than the heat sink. The isotherms are parallel 
for depths below 50 pm. 
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Temperature profile through AIGaN, GaN, AIN 
SiC (V d =20 V) 



Figure 85. — Temperature profile though the upper portion of the HEMT. The AIGaN and AIN layers are on the left and 
right sides in the right panel. 


The final group of figures investigates the “singularity” behavior just under the gate edge. Figure 86 
displays a very thin portion of the heterojunction; the contours are those of total current density. The spike 
region is extremely devoid of current. Notice it extends below the transition layers which means the 
“channel” no longer exists along the heterojunction. In other words the carriers are in a channel to the left, 
then they dip down into the bulk GaN (go around the notch) then return to the heterojunction channel on 
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the right. Figure 87 shows this effect from a wider vantage point. The horizontal outline (designated 1 in 
the left panel) demonstrates the reduction of current at the gate edge plane (x = 52 pm) as seen in the right 
panel. The value at the notch is 4000 A/cm 2 , which is three orders of magnitude below that on either end 
of the outline. Figure 88 shows the conduction band edge along this plane, and observe the notch has been 
removed by the electric field due to the gate and drain potentials. Figure 89 is a closer view of the band 
edge which demonstrates the effective dissolution of the channel in this region. Figure 90 illuminates the 
behavior of the magnitude of the electric field. The value at the notch is 4.6x 10 4 V/cm. The panel on the 
left is a horizontal cut at the plane y = 0.0502 and 50.4 < x < 53.4. While that on the right is a vertical one 
withx = 51.7, 0.0494 <y< 0.0535. The step in the field is due to the sheet charge placed at the 
heterojunction. 

Finally figure 91 demonstrates the singularity in another fashion. The two vertical cutlines are given 
in the left panel. The middle panel shows the carrier density at the heterojunction to be about 2x 10 14 /cm 3 
at the singularity. Flowever, at the midpoint of the access region it is about 4xl0 19 /cm\ which is typical of 
that in the channel in most devices. 
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Detail of the GaN channel-to-bulk transition region 
SiC {Vd = 20 V) 
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Figure 86. — Close view of the AIGaN/GaN heterojunction plane at 
0.05 pm and the thin transition layers separating the “channel” 
properties from the bulk GaN properties. 
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Electric field under gate edge 
SiC (Vd= 20 V) 




Figure 90. — Variation of the magnitude of the electric field in the singularity region. The left and right panels are 
horizontal and vertical cuts respectively. 
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Carrier Density under Gate Edge 
SiC (Vd= 20 V) 



Figure 91. — The carrier density at the gate edge and midway the gate-drain access region. Notice the carriers in 
the singularity region are about 5 orders of magnitude below the normal value contained in the conduction 
band notch of a typical heterojunction. 


X. Multiconductor Results 

This section covers the multi-finger comparisons with Kuball’s data (ref. 17). The measured 
temperatures and schematic of the device are given in figure 92. The dotted curves that blend with the 
measured solid peaks were derived from a thermal simulator. The device is designated as the 8x250 pm 2 
unit with 25 pm spacing between gates. Table 17 displays the measured temperatures on all eight gates. 
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TABLE 17.— TEMPERATURES ON THE GATES OF THE 8x250 mm 2 DEVICE. 

THE GATES ARE 0.8 pm LONG, 250 pm WIDE, THE SOURCE-DRAIN 
SPACING IS 5 pm, AND THE SUBSTRATE IS SiC 


Gate no. 

1 

2 

3 

4 

5 

6 

7 

8 

T, °C 

207 

215 

233 

224 

234 

215 

218 

206 


Close scrutiny reveals an interesting pattern; the two gates in the center (numbers 4 and 5) differ by 
10°. Yet gates 3 and 5 are at the same temperature. This either demonstrates the uncertainty in the 
measurements or a nonuniformity in the device. At any rate, we can say the outermost gates are at the 
same temperature and those in the center are about 20° hotter. 

In table 18 we present the ratios of temperatures; the second column gives that for the measured 
values; while columns 3 and 4 display the predicted values using either the distributed or partial 
inductance approaches. 


TABLE 18.— RATIOS OF GATE TEMPERATURES; MEASURED, AND PREDICTED 
BY THE TWO INDUCTANCE METHODS (DISTRIBUTED AND PARTIAL) 


Ratio 

Measured 

Distributed 

inductance 

Partial 

inductance 

T 2 /Ti 

1.04 

1.09 

1.03 

T->JT\ 

1.13 

1.14 

1.04 

m 

1.08 

1.16 

1.05 

Js/r, 

1.13 

1.16 

1.05 

T(JT\ 

1.13 

1.14 

1.04 

T-jITi 

1.05 

1.09 

1.03 

zyr, 

1 

1 

1 



— -i h*— 0.8 p 

] □ 1 = 

H* — 5 p— H 



Figure 92. — Kuball’s measured temperatures (°C) on an eight gate device on SiC. 
The perspective shows the gates numbered from left to right, along with the 
dimensions of the chip. The gate length and S-D separation are 0.8 and 5 pm 
respectively. 
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Figure 93. — Display of the ratios of measured temperatures 
and the corresponding ratios for the distributed and 
partial inductance approximations. 7(5)/7(1) is the ratio 
of the temperatures on gates 5 and 1 . 


Inspection of the first column shows that a nearly symmetrical pattern exists, and some discrepancy 
exists for gate 4 (it is too low). Comparison of all three columns shows that the distributed approach gives 
a stronger variation, but the scatter in the data tends to make the results fall between the distributed and 
partial inductance approaches. Thus either approach will predict the actual ratio to within 1 to 10 percent. 

Recall the prescription given in section IV for multiple fingers assumes the Rn matrix element is 
known, and it is assumed to be the thermal resistance of a single gate on the given substrate. The 
remaining elements are then determined from scaling by the factors of the mutual inductances between 
the gates. Therefore, to have any chance at good correlation, the value of Rn must be accurate. We arrive 
at Rn in the following manner. First we use the results from the simulation. At 20 V the simulation found 
the drain current to be 91 mA, so the dissipated power was 1.82 W. Table 19 provides a summary of the 
temperatures determined for a single stripe on SiC. 


TABLE 19.— TEMPERATURES AT THE INTERFACES FOR THE CASE OF 
A SINGLE STRIPE ON SiC. THE DISSIPATED POWER WAS 1.82 W 


Location 

Base 

SiC, 

top 

AIN 

GaN, 

top 

AlGaN 

T, °C 

32 

123 

126 

148 

148.5 

T, K 

305 

396 

399 

421 

421.5 


The actual temperature on the top surface of the SiC is 123 °C. The corresponding linearized 
temperature 0=96 °C. In the simulation we used k(7) = 5.2(300/7) 2 4 . Thus the thermal resistance, 

Rn, = (96 - 32)/l .82 = 35.2 K/W. The gate length is 1 pm, and if we use the footprint that results from 
assuming 45° spreading in the GaN, we have the value 3.4 pm for the source on the top surface of the 
SiC. Table 20 displays the thermal resistances versus the gate length 2l x . 


TABLE 20.— THE CHANGE IN THERMAL RESISTANCE WITH 
GATE LENGTH FOR A SINGLE STRIPE ON SiC 


2 1„ pm 

3.8 

3.4 

3.2 

0.5 

0.25 

0.125 

R th , K/W 

18.8 

19.3 

19.6 

28.5 

31.9 

35.3 


Notice that the last entry gives the value determined from the simulation. Similarly, the actual 
temperature on the GaN is 148 °C, with corresponding linearized value being 0 = 147 °C. The assumed 
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variation of kappa was k( 2)= 1.6(300/7)' 4 . ThenTG, = (147 - 123)/1.82 = 13.2 K/W. By sweeping 21 x we 
obtain R th = 13.6 at 21 x = 3 pm. This is exactly the gate-drain separation. 

From the above calculations we may draw the following conclusions. The best way to determine R\\ 
is to use the temperatures found from the simulation; thereby eliminating the step of attempting to 
estimate the size of the effective heat source. The procedure is to start with the temperatures developed by 
simulation, then use the Kirchhoff transform to find the corresponding linearized temperatures. Also from 
the simulations one may determine the heat flux into the sink by using the equations and method outlined 
in the companion memo. Then calculate Rn from its definition; change in linearized temperature divided 
by heat flux. The tables above have shown that the “effective” heat source is some value between a 
fraction of the gate length and the source-drain separation. Flowever, the final calculated temperature is 
too sensitive to the chosen value to be reliable. Once Rn is determined from the 2-D simulation, then the 
multiconductor calculation will yield very good values for gate temperatures for multi-gate devices. 

XI. Discussion and Conclusions 

The modeling of AlGaN/GaN FIEMTs is not easy as there are many gray areas that have not been 
resolved at this time. Due to this situation, the present work necessarily had many arbitrary assumptions. 
We discuss many of the problematic points that can/should be further investigated to improve the 
modeling of this promising technology. 


Problematic Issues 

1 . The AlGaN/GaN F1EMT is a semiconductor device that couples properties of heterojunctions, 
mechanical strain, piezoelectric charge, thermal and electric fields. A simulator needs to be developed 
that self-consistently solves for all of the above mentioned constraints. At the present time (August 2008) 
none exist. 

2. We observed the sensitivity of temperature with k( 7), and this must be known accurately for all 
layers. Apparently the growth and processing steps significantly alter this parameter, so this must be 
quantified for future work. The representation for k(7) of SiC was expressed as the reciprocal cubic as 
stated in section VII. The (k,7) pairs chosen are given on page 35; the temperatures selected were 350, 
400, and 450 K. When extrapolated back to 300 K, the formula yields 5.24 W/cm-K for k. The SiC for 
Kuball’s devices were from Cree (the only source they mentioned), and measurements from Cree indicate 
k is near 4.85 at 300 K for 4H-SiC (SI). Thus our best fit used a value (5.24) that was 8 percent high. 
There are several reasons for this; (a) we assumed the baseplate remained at 300 K, whereas it actually 
drifted up in temperature, and there is no way to allow it to drift upwards during a simulation run. Thus to 
get correlation with their measured points, we had to use a slightly higher value for k. (b) The device had 
some heat escaping from the sides and metal pads on the top; these heat paths will lower the gate 
temperatures from those under the assumption that all heat must leave via the baseplate, (c) We could 
have chosen the 300 K value for k in the curve-fitting, rather than allowing the extrapolation back to 

300 K to perhaps appear out of bounds. In the literature one often finds 3.3 to 3.5 as values for k, but that 
is for partially conducting material (61 I-SiC) where the doping is near l.OxlO 17 cm” 3 , not SI (semi- 
insulating). (d) Part of the 8 percent is due to the 2-D and 3-D effective differences. Therefore we are 
comfortable with the curve-fit chosen. 

3. The papers by Kuball’s group are the only ones in the open literature with enough data to permit 
validation of modeling with experiments; more are needed. We could not make an unambiguous 
comparison as the baseplate temperature was not fixed. We chose to keep the gate in the middle of our 
simulation space, although the experimental gate was not centered. We assumed this difference to be 
inconsequential. 

4. A significant problem is the source of the carriers in the notch of the conduction band. At the 
present time there are several theories that address this issue, but none have been justified. Without 
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knowing the true fixed charge density, Poisson’s equation cannot be solved; which is the heart of the 
semiconductor simulator. We used a supply layer in the AlGaN layer, and the mobile charges inserted 
along with the sheet charge placed at the heterojunction. ATLAS automatically inserts them to keep 
charge neutrality. Many practical devices are unintentionally doped (UID), so the source of the carriers is 
u nk nown. The source for (some of the) carriers has been speculated to be donors on the upper surface of 
the AlGaN. Most (if not all) devices require a passivation layer on the top surface, which chemically 
interacts with the donors (neutralizes them?). 

5. The device has many traps, and they are apparently distributed in all layers and interfaces. The 
interplay between traps, surface donors, and free carriers in the 2DEG channel cause considerable 
variation while in operation. The phenomena of current collapse, rf slump, knee walk-out, and other non- 
ideal effects are assumed to be related to the traps, u nk nown donors, and other effects. Current collapse is 
assumed to be due to traps in the GaN, whereas rf slump is due to surface traps. 

6. We used a very simple model for the 2DEG channel; just a layer of GaN with appropriate mobility. 
We performed a calculation using the Poisson- Schroedinger solver, and determined the I-V curve for 0 V 
on the gate. The curve was close to what we obtained using the simple model, so it is justified. 

7. The density of carriers in the notch, n s , apparently change with many design parameters. The 
interaction is complex as gate length changes (all other parameters held constant) strongly alter the 
saturation drain current (experimental finding). Since the source-drain separation is constant, apparently 
the amount of exposed surface area between them affects the current. One speculates this may be due to 
the change in the number of donors or traps (not covered, and hence neutralized by gate metal) with 
change in available AlGaN surface area. To account for this and other effects (for example the change in 
strain with applied voltage and the corresponding change in polarization charge density, and thus change 
in 2DEG charge density), the following prescription could be used. Write the 2DEG carrier density as 
(see refs. 80 and 8 1 ) n s =j{V gs ,T), where to first order, n s varies linearly with V gs . Since the simulator 
charge sheet has a strong affect on the 2DEG density, one has the freedom to vary it to change n s for 
different gate voltages. This could act along with changes in the supply layer to generate the appropriate 
2DEG density. This would enable better correlation with SiC than we were able to achieve. The device on 
the sapphire substrate was larger and hotter than that on SiC, and the same value for a„ was sufficient for 
all gate values. However, the SiC device, being smaller, was more sensitive to gate voltage. Another 
possible option could be the fact that the sapphire unit, being hotter, had most, if not all traps emptied; 
whereas the cooler SiC had traps remaining, which caused some of the droop. By adjusting the sheet 
density, we are indirectly taking into account the reduction in channel current due to charge now held in 
traps. However, we are not directly including this trapped charge in Poisson’s equation. 

8. Many semiconductor properties were not included; presence of holes (some simulations in the 
literature, due to their arbitrary assumptions, predict a substantial number), no tunneling between 2-D and 
3-D (normal) electrons in and out of the well, and no traps. The presence of traps on the surface may pin 
the Fermi level; this was not considered. Also some of the droop in the I-V could be due to trapping. No 
hot electrons were included since the D-D model was used. However, some experts in the field think hot 
electrons don’t exist in the devices. This is based on experimental data. 

9. Basically we tweaked about 14 basic material parameters to achieve satisfactory agreement with 
the simulator and measurements. This is needed, as many of them are uncertain by up to 30 percent, as 
well as the shortfalls of the program; i.e., no strain, thus the piezoelectric effects are masked in just a 
uniform sheet charge, and a 2-D, not a 3-D self-consistent solver. The supply layer in the AlGaN serves a 
critical need in a simulation. In actual devices the carriers that populate the channel may come from either 
a supply layer or unintentional dopants in the AlGaN layer. In either case, dopants must be present. Three 
publications (refs. 73 to 75) have argued that the electrons in the channel cannot come from the bulk, and 
therefore must come from intentional or unintentional dopants in the AlGaN. The dopant density is a 
dominant parameter, as it serves to adequately supply the notch with sufficient carriers. Relatively small 
changes in the supply layer density produce enough carriers to cause the I-V curves to rise steeply. In the 
UID devices reported, there has been speculation that the dopants are on the upper surface of the AlGaN. 
Then the details of the passivation layer become important; as it may restrict the dopant density and thus 
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the available drain current. The passivation layers stabilize the drain current with time and temperature; 
and we speculate that the mechanism is essentially the changing of the dopant density activation level. 

The number of available carriers adjusts the shape and depth of the notch in the conduction band edge, 
and therefore the 1-V curves significantly. In one simulation, a negative sheet charge was placed on the 
top AlGaN surface, and the drain current was reduced by 40 percent of that without the sheet. The sheet 
was assumed to model the negative spontaneous polarization charge layer that must exist at the surface. 
However, apparently in actual devices this charge is neutralized with positive charge, due to dopants or 
other u nk nown sources. The interaction of the surface/passivation layer interface and the channel 
parameters is a topic for future investigation with ATLAS. 

Since experiments have shown that irreparable damage occurs when the channel approaches about 
230 °C (or 600 °C, depending on reference), our model can be used to determine the maximum DC power 
that can be dissipated safely for various designs. Also the singularities about the gate edge can be studied 
for both excess temperature and field (to eliminate the possibility of avalanche) simultaneously. Then the 
maximum drain voltage (and hence maximum RF power theoretically obtainable) may be estimated. We 
attempted to use the Barnes model (a choice for negative differential mobility) in the SiC cases, as we had 
difficulty reproducing the large experimental droop. We did not get satisfactory results; perhaps 
instability problems were occurring (a warning of such was stated in the users’ manual). 

Figure 94 gives the entire energy band (from gate to heatsink) for both substrates. The values of 
energy in eV at selected points are referenced to the Fermi level. The slope in the SiC region is due to the 
Schottky barrier height being applied to the bottom of the SiC (since in contact with a metal), as SiC is 
recognized as a semiconductor in ATLAS. This contact potential should not alter our results; however its 
effects may be seen in some plots, wherein the potentials in the drain area are a few volts above 20 V, 
even though the drain is biased to 20 V. The sapphire is treated as a dielectric and thus has no contact 
potential. 
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Figure 94. — Energy band profile through entire device for SiC and sapphire substrates. 
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The conduction band edge shape in the vicinity of the heterojunction is one of the major components 
that determine the characteristics of a HEMT. The step height and shape of the notch determines many of 
the 2DEG properties. The step is due to the difference in bandgaps in the materials on either side of the 
junction. The change in lattice parameter is the fundamental cause for the difference. At the material 
interface the two space lattices must transition over a short distance, and the details of the transition are 
not well known. The following calculation determines the effective dipole needed to model the step if the 
step were due to classical static charge layers. From (ref. 76) we showed the step, or more accurately, the 
linear drop across two infinite sheet charges that form a dipole to be 

A(j) = ad/s 

Where a is the magnitude of the sheet charge density on one sheet, d is the separation between them, and 
s is the dielectric constant. For our case, the step is 0.45 V, 8 = 9.5 s 0 , and <7= 1.92x10 8 cm, which yields 
a = 1 .23 x 1 0 14 electrons/cm 2 . The separation was approximately the distance between three mesh points. 
Observe that this is about two orders of magnitude larger than the sheet charge used at the heterojunction 
to model the piezoelectric fields. The step, being of quantum mechanical origin, apparently has no 
effective charge associated with it, so only the real positive sheet charge and the mobile carriers in the 
notch create electric fields in the region of the heterojunction. Thus our model implies a “clean” 
heterojunction with no free or bound charge due to the lattice change. The only bound charge is that 
which creates the piezoelectric field, and it exists in all regions in the vicinity of the heterojunction where 
the crystal is strained. We have shown that this distributed effect may be adequately modeled with a 
single sheet charge at the heterojunction. The electric potential field associated with the bound charges 
due to strain, along with other terms mentioned in (ref. 76) constitute the “other” terms that constitute the 
conduction band edge energy profile. 

Figures 95 through 99 display the device current at several vertical planes. The vertical outlines are at 
the center of the S-G access region, the gate edge, and the center of the G-D access region. The total, 
channel, AlGaN and bulk layer currents are displayed. The supply layer doping value is 2x 10 18 cm 3 . 

Figures 100 through 103 are similar plots, but now the doping is 3xl0 18 cm' 3 , and the gate is at -2 for 
all cases. 


Division of current ( V s - 0, Vd = 4) 
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Figure 95. — The division of current at three vertical outlines. Observe the 
total value should be absolutely flat, but the small variation is due to 
deviations in the numerical integration of the current density due to 
non-uniform mesh densities. 
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Division of current (V g = 0, Vd = 8) 



Figure 96. — Current division at 8 V. 
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Division of current ( V g = 0, Vd = 12) 



Figure 97. — Division of current at 12 V. 
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Current (mA) Current (mA) 


Division of current (V g = 0, Vd = 20) 



Figure 98. — Division of current at 20 V. 
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Figure 99. — Division of current for V g = -2V and Vd = 20 V. Notice all of the current is 
in the bulk at the gate edge (52 gm). 
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Figure 101. — Division for drain at 8 V. Notice as the drain increases, the channel 
becomes depleted at the gate edge. 
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Figure 102. — Similar division at 12 V. 
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Division of current ( V g = -2, Vd = 20) 
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Figure 103. — Current division for SiC device at bias point (-2,20). The channel is 
completely empty of carriers at the gate edge. 
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The above sequence of plots shows that the 2DEG channel is devoid of carriers over a small region 
beyond the gate in the saturation region. Thus the conductivities of the bulk GaN and AlGaN layers are 
important in device design. 

We have used Masana’s equation set as the basis of our analytical efforts, since it appeared to be the 
most flexible. However, by referring to figure 1 04 we observe the large variation of thermal resistance 
with the “effective gate length 21”, for most of the closed- form expressions published in the literature. 
The legend contains the works by Meeks (ref. 77), Cooke (ref. 78), and Fukui (ref. 79), along with others 
that we have already discussed. While all show a decrease in R th with increasing 2 l„ the spread is too 
large if one is attempting to develop accurate results. For a 0.5 pm gate (heat source) we observe a spread 
of about 3:1 from Royet’s to Darwish’s equations. Notice also, that Darwish’s results are not applicable 
for heat sources greater than 1.5 pm. With this new insight, we decided that our first attempts to find 
“effective” gate values, as presented in section V, with many reservations, was not the best way to 
proceed. 

We have concluded that one must abandon the “effective gate length 21” approach, and determine 
instead, R,i, directly from the simulations. That is we determine R th from its definition; change in 
temperature between the planes of interest, divided by the power dissipation. We found this to be 
necessary, as the heat generation region is not adequately modeled by a uniform heat source strip with the 
dimensions equal to the actual gate, or some constant times the gate length, or some value near the 
source-drain spacing. The simulations show that there are in general 4 spikes in the heat generation 
function; which are located at the metallization edges. Their relative intensities change with bias point, 
which makes determining an effective length for all reasonable bias points nearly impossible. In addition 
there is the distributed portion that also changes with current path through the device. The details of the 
determination of R th via the simulation will be presented in another memorandum. Also presented there in 
the appendix is the determination of the heat flux q between the ends of a 1-D flow regime with known 
temperatures at either end, and the thermal conductivity described by an easily integratable function of 
temperature. This is necessary, as q is not immediately available from ATFAS. We found that q may be 
only about 60 percent (with a sapphire substrate) of the I-V product of the bias point. This gives an 
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Figure 104. — The thermal resistance versus the effective gate length 
2/ x , for seven of the basic references that calculate it in closed-form 
equations. This chart should be nearly exhaustive, as the literature 
was probed to great depth. 
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indication of the amount of heat that leaves via the sides, and metal pads on the surface. Notice also 
that both components in equation (1) R th and Pd iss , (where it is assumed that P diss = q ), are not accurately 
found from formulas or the 1-V product of the bias point! This is a significant result obtained in this 
study. 

Accurate temperature predictions are critical, as lifetime estimates can change by large values as the 
temperature changes by 10°. While we were able to get agreement, there are areas that leave the door 
open for speculation. First of all the actual devices had 3-D heat spreading as opposed to the 2-D 
simulation. For example, the top surface of the device had dimensions 297x179 pm 2 . The 200 pm wide 
gate was 70 pm from the edge. Flowever, in the simulations, the areas were 120x200 (sapphire) and 
1 14x200 (SiC) pm 2 , and the gate was symmetrically located in the center. This was due to the upper limit 
of mesh points available in ATLAS. In the simulation the heat can only flow perpendicular to the long 
dimension of the gate electrode, and not “along” the gate. The short sides of the simulations were dictated 
by the number of available mesh points and the maximum separation between them needed for accuracy. 
Thus the ratios of the true surface area to the simulation areas were 2.2 (sapphire) and 2.3 (SiC). Thus the 
actual devices had more area upon which to dissipate heat. Thus we realize our values for k are too large; 
however, once the effective value is determined (to compensate for the discrepancy in areas and the lack 
of spreading along the gate width), then the simulations tracked the measured values very well. The 
anticipated future work includes the following tasks. The incoiporation of traps will be the next major 
hurdle in the simulation effort. The behavior with traps is a complex and not well understood 
phenomenon. The effects of trapped charge alter both the DC and rf characteristics in complex ways. We 
have studied and summarized the literature concerning physical locations, excitation levels, surface and 
bulk densities, and fill-empty times for all traps that have been measured. Once the trap incoiporation is 
completed, we will be able to accurately determine the benefits, drawbacks, etc. of alternate gate shapes, 
single or double recess gate designs, fieldplates, and passivation layers. Another exercise is to compare 
results of ATLAS 2-D (with just the thermal simulator) and ANSYS 3-D to quantify the differences in 
temperature. One can use both uniform stripes as well as the more realistic four spikes and a smooth 
portion of g(x,y). Favorable agreement could resolve the magnitude of the error in the 2-D case. 

XII. Summary 

This memo has described the development of a baseline model for simulation of AlGaN/GaN HEMTs 
with the ATLAS package from Silvaco, Inc. We have sought to understand the basic physical principles 
underlying the device, and apply the available tools in ATLAS to allow reasonable simulation. We used 
the enhanced drift-diffusion model, as experiments have indicated that hot electrons are not plentiful in 
the devices, and the hydro-dynamic model is not necessary. A single sheet charge at the heterojunction is 
sufficient to model most of the stress/strain, and perhaps trap effects. The channel supporting the 2DEG 
was modeled with a uniform slice of GaN with appropriate mobility parameters. A supply layer is needed 
in the buffer to guarantee sufficient carriers in the device. Apparently the carriers in the 2DEG channel 
are insufficient to carry all of the measured drain current, (this fact has been stated in the literature). We 
have demonstrated many of the sensitivities between various parameters in the system. The 1-V and 
thermal characteristics of single gate devices on sapphire and SiC substrates were developed. The self- 
heating included the dependence of the thermal conductivities of all materials. The 2-D simulations for 
single finger devices agree very well with published data in both the I-V characteristics and channel 
temperatures. Therefore we think the basic model and parameter set is adequate for predicting I-V and 
temperature values in all part of the transistor as device parameters and geometry are changed. 

We have presented and used data from published measurements, simulations, and theoretical studies, 
to develop the channel temperature estimation procedure. We strove for self-consistency, as the 
dispersion in parameter values is very large. This dispersion is to be expected due to the various growth 
methods, processes, and other procedures. Our simulations have shown that the current in the channel is 
the dominant portion in the linear regime of the I-V curve, and then in saturation, the channel is severed 
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by the large field at the gate edge, and the current is then carried in the bulk GaN and AlGaN. We assume 
the step height in the conduction band at the heterojunction includes any variations due to dipoles, excess 
charges, etc. While the details of the notch in the band edge are important, the contact resistivity for the 
source and drain contacts, and the doping in the AlGaN supply layer are also very critical parameters. For 
our software version, we were at the limit of the number of mesh points; which necessitated the 
simulation of just a single gate device. Since the simulator was restricted to 2-D calculations, the 
analytical portions that include 3-D effects was needed to complete the study. We have determined that 
the major difficulty in matching analytical predictions with measurements resides in choosing the 
appropriate size of the heat source in the established closed form equations found in the literature. The 
effective gate length must adjust as the heat generation g(x,y,z ) changes with bias. The change in thermal 
conductivity with temperature must be included in any simulation or analytical calculation or significant 
error will result. The channel temperatures are dominated by the substrate, since the GaN region is very 
thin compared to the substrates. 

The electric fields and electron velocities were very large at the contacts; with the source being the 
larger. The fields were up to 7 MV/cm, which are slightly above the accepted avalanche values of 3 to 
4 MV/cm. The carrier velocities were as high as 4x 1 0 ? cm/sec, even with the designated saturation 
velocity set at 2x 10 7 cm/sec. These minor effects were assumed unimportant, and did not detract from the 
overall favorable outcomes. We found the channel carries about Vi of the total drain current; with the 
remainder nearly split between the AlGaN and bulk regions. In a small region about the gate edge on the 
drain side, the channel is broken for voltages above 12 V. It recovers at somewhere between the gate edge 
and the middle of the G-D access region. We have demonstrated that the analogy between the heat 
conduction and Poisson’s equation may be used for developing the thermal resistance matrix. The effects 
of finite width gates can be incorporated using the partial inductances between the gates. It was 
demonstrated that the partial inductance procedure gives a less severe change in temperature between 
gates; which seems physically correct, as heat is allowed to escape along the gates, which should bring 
them all to lower as well as closer temperatures. 

The major predictive nature of is work is as follows. Suppose one wants to construct a multi-finger 
power device. First construct a single finger device and measure the I-V curves. Use the simulator to 
duplicate those curves. Once this step is completed, the temperatures at different regions in the device are 
known. Determine the thermal resistance for the single finger from the results of the simulation; 
temperature change divided by the heat flux. Then R n is known, and the remaining analytical steps have 
been discussed herein. Then the temperature on each finger of the predicted array is calculated. This 
enables design changes with minimal effort, and circumvents the very difficult task of accurately 
measuring the temperature. The accuracy of the technique should be within 5 to 7 percent. 
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